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SECTION  I 


INTRODUCTION 


DYNAMIC  MATERIAL  BEHAVIOR 

The  factors  common  to  any  static,  or  dynamic,  stress 
analysis  problem  consist  of  the  following:  the  specimen 
geometry,  the  loading  applied  at  the  boundary,  and  the  material 
of  the  specimen.  These  three  factors  will  interact  to  produce 
the  stress  level  inside  the  body.  The  response  of  the  body  to 
these  stresses  is  important  in  most  engineering  endeavors. 
Frequently,  a  rational  engineering  design  requires  the  ability 
to  predict  the  stresses  and  material  response  in  each  component. 

The  most  common  approach  to  investigating  the  response  of  a 
material  is  to  fix  two  of  the  three  factors  governing  the  stress 
problem.  As  an  example,  in  a  uniaxial  tensile  test  the  specimen 
geometry  and  loading  have  been  standardized.  This  provides  a 
uniform  basis  to  compare  the  third  factor,  the  specimen  material. 
This  test  is  based  upon  a  straightforward  method  of  calculating 
the  stress  level  in  the  specimen. 

There  are  several  limitations  to  this  type  of  materials 
investigation.  First,  the  response  of  the  material  to  the 
complex  loadings  experienced  in  service  may  not  be  accurately 
represented  by  such  a  simple  test.  Second,  for  some  applications 
the  experimental  apparatus  cannot  provide  the  necessary  loading 
requirements  seen  in  service,  e.g.,  high  strain  rates. 

To  overcome  these  restrictions  the  designer  can  take  either 
of  the  two  following  approaches:  use  a  large  factor  of  safety, 
or  proceed  on  a  need  to  know  basis.  The  former,  though  widely 
used,  will  not  be  discussed.  The  second  approach  usually  begins 
by  redefining  the  experimental  technique  in  terms  of  the  specific 
problem  at  hand.  This  approach  is  used  frequently  for 
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investigating  materials  response  to  rapidly  applied  loads.  If 
the  loading  rate  and  magnitude  are  sufficiently  high,  the 
material  response  is  called  dynamic. 

Dynamic  material  behavior  has  been  characterized  by  the 
presence  of  inertial  effects  and  wave  propagation  which  affect 
the  stress  distribution  inside  the  specimen.  If  an  impulsive 
load,  large  enough  to  cause  permanent  deformation  is  applied  to 
the  boundary  of  a  specimen,  the  stress  in  the  region  nearest  the 
load  will  be  significantly  higher  than  in  any  other  portion  of 
the  body.  The  deformation  which  occurs  in  the  specimen  can  be 
modeled  by  its  wave-like  motion  Through  the  material.  If  the 
deformation  takes  place  rapidly,  the  particle  being  displaced 
will  have  some  inertial  energy.  If  the  inertial  energy  is  large 
enough,  it  can  have  a  significant  effect  on  the  final 
configuration  of  the  specimen. 

Certain  aspects  of  dynamic  behavior  have  been  known  since  the 
19th  century.  British  investigators  showed  that  an  iron  wire 
could  resist  permanent  deformation  under  large  loads  for  short 
periods  of  time.  This  test  proved  that  a  relationship  existed 
between  the  yield  stress  in  a  material  and  the  rate  at  which  the 
load  was  applied. 

In  the  1940's,  G.I.  Taylor,  Reference  1,  met  with  some 
success  at  charact  izing  this  behavior.  He  proposed  an 
experiment  to  measure  what  was  then  called  the  dynamic  yield 
strength . 

The  experiment  proposed  by  Taylor  has  become  a  standard  test 
in  laboratories  that  study  the  behavior  of  materials  at  high 
rates  of  deformation.  The  Taylor  test  consists  of  impacting  a 
plane-ended  cylindrical  projectile  against  a  relatively  rigid, 
massive  anvil.  What  should  come  out  of  the  Taylor  test  is  the 
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yield  stress  level  for  a  material  that  is  rapidly,  or 
impulsively,  loaded.  The  Taylor  test  was  designed  to  standardize 
the  specimen  geometry  and  the  loading  pattern  applied  to  the 
specimen  boundary.  As  previously  mentioned,  the  last  component 
of  the  internal  stress  problem  would  be  the  material  under 
investigation . 

One  use  for  the  information  generated  by  the  Taylor  test  is 
in  the  development  of  armor  and  armor  penetrators.  Models  of  the 
interaction  between  a  target  and  a  penetrator,  References  2  and 
3,  require  that  the  materials  be  characterized  by  their  dynamic 
strength  values.  Another  use  for  the  Taylor  test  is  as  an 
accuracy  check  of  two-dimensional  computer  models  of  deformation 
behavior,  References  4  and  5. 

The  object  of  this  report  is  fourfold.  First,  it  is  to 
examine  aspects  of  both  two-dimensional  and  one-dimensional 
modeling  used  with  Taylor  testing.  Second,  it  is  to  describe  in 
detail  a  recently  constructed  Taylor  test  apparatus.  Third,  it 
is  to  provide  an  analysis  of  data  obtained  experimentally  and  . 
used  in  the  one-dimensional  models.  Fourth,  it  is  to  report  on 
observations  made  on  the  microstructure  found  in  test  specimens. 
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SECTION  II 


ONE-DIMENSIONAL  MODELS 


1 .  BACKGROUND 

One-dimensional  models  of  the  Taylor  experiment  are  used  to 
calculate  the  yield  stress  level  of  a  material  from  post-test 
measurements  of  specimen  deformation.  Historically,  all 
one-dimensional  models  are  based  on  the  analysis  developed  by 
Taylor,  Reference  1.  Over  the  years  various  investigators  have 
proposed  modifications  to  his  analysis  by  changing  the  basic 
equations  or  using  different  types  of  material  constitutive 
relationships.  Therefore,  the  discussion  of  one-dimensional 
models  should  begin  with  a  development  of  Taylor's  original 
analysis  of  the  problem. 

At  impact,  a  wave  of  compressive  stress  will  be  generated  at 
the  anvil  face.  If  the  velocity  of  the  projectile  is 
sufficiently  high,  the  stress  wave  will  separate  into  two 
components.  The  first,  or  leading,  component  is  an  elastic 
compressive  wave,  moving  through  the  material  at  the  speed  of 
sound.  The  amplitude  of  the  stress  level  behind  the  compressive 
wave  front  is  below  the  yield  strength  of  the  material.  The 
second  component,  a  plastic  compressive  wave,  will  follow  the 
elastic  compressive  wave  at  a  greatly  reduced  velocity.  At  the 
plastic  front  the  stress  level  exceeds  the  yield  strength  of  the 
material.  The  high  compressive  stress  causes  severe  deformation 
to  occur  in  the  form  of  radial  motion  outward  away  from  the 
specimen  axis,  accompanied  by  axial  shortening  of  the  specimen. 

As  the  event  proceeds,  the  elastic  compressive  wave  will 
arrive  at  the  free-end  of  the  specimen,  where  it  is  reflected  as 
a  tensile  wave  of  equal  magnitude.  The  tensile  wave  will  move 
through  the  specimen  until  it  encounters  the  plastic  compressive 
wave  front,  located  within  the  specimen,  Figure  1,  Reference  6. 
The  motion  of  the  elastic  wave  and  the  interaction  with  the 
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plastic  compressive  wave  will  have  two  important  effects.  First, 
the  velocity  of  the  plastically  undeformed  portion  of  the  rod 
will  be  reduced  as  the  elastic  wave  moves  through  the  material. 
Second,  the  reflected  tensile  wave  will  superimpose  with  the 
compressive  plastic  wave  to  reduce  the  overall  stress  at  the 
plastic  wave  front.  After  repeated  occurrences,  the  motion  of 
and  stress  within  the  specimen  will  both  be  reduced  to  zero. 


PLASTIC  WAVE 
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Figure  1.  Plastic  and  Elastic  Wave  Motion  in  the  Specimen 


2.  TAYLOR'S  ANALYSIS 

To  construct  a  model  of  the  impact  event,  Taylor  makes  three 
assumptions  in  his  analysis:  the  material  stress-strain 
relationship  is  rigid,  plastic;  radial  inertia  effects  can  be 
neglected;  and,  a  condition  of  uniaxial  stress  exists  across  the 
elastic/plastic  interface.  The  relative  effects  of  these 
assumptions  have  stirred  numerous  debates  and  papers  regarding 
the  validity  of  his  analysis.  The  simplicity  of  the  experimental 
technique  and  subsequent  reduction  of  data  are  incentives  to 
accept  these  assumptions.  It  must,  however,  be  kept  in  mind  by 
the  user  of  such  data  the  level  of  approximation  that  was  used  in 
the  construction  of  the  analysis. 

Taylor's  analysis  relates  the  altered  geometry  of  the 
specimen  after  impact  to  the  dynamic  yield  strength  of  the 
material.  In  this  way,  Taylor  could  extract  the  crucial  dynamic 
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strength  from  only  two  postmortem  measurements  of  the  deformed 
specimen.  Taylor  formulates  his  analysis  through  equations  that 
relate  various  kinematic  parameters  during  the  impact  event,  such 
as  the  time  required  for  an  elastic  wave  to  travel  down  the  rigid 
portion  of  the  specimen  and  back  to  the  plastic  wave  front,  the 
incremental  change  in  the  position  of  the  plastic  wave  front,  the 
foreshortening  of  the  rigid  portion  of  the  rod,  and  the 
incremental  change  in  the  velocity  of  the  undeformed  portion  of 
the  rod.  By  eliminating  the  speed  of  sound  in  the  material,  he 
generates  a  set  of  differential  equations.  These  differential 
equations  define  the  velocity  of  the  plastic  wave,  the  rate  of 
foreshortening  of  the  undeformed  portion  of  the  rod  and  its 
deceleration . 

To  begin  the  development  of  the  analysis  by  Taylor,  it  is 
first  necessary  to  define  the  nomenclature  to  be  used,  Figure  2. 
Let  L  represent  the  original  length  of  the  specimen,  and  S,  the 
time  dependent  displacement  of  the  undeformed  portion  of  the  rod 
relative  to  the  initial  configuration.  At  some  time  after 
impact,  X  represents  the  extent  of  the  plastic  zone  relative  to 
the  original  configuration.  The  position  of  the  plastic  front 
is  h,  measured  relative  to  the  anvil  face.  The  current  length 
of  the  undeformed  portion  is  given  by  1 .  A  relationship  exists 
between  the  time  dependent  quantities,  S,  1,  h,  and  the  original 
length 

L  =  S  +  i  +  h  (1) 

Differentiating  Equation  (1)  with  respect  to  time  gives 

•  •  • 

o  =  s  +  a  +  h  (2) 

or 

•  •  • 

h  =  -(S  +  a)  (3) 

But  S  is  simply  the  velocity,  v,  of  the  back  end  of  the  specimen, 
and,  h  is  renamed,  \,  the  Eulerian  plastic  wave  speed,  to  give 

\  =  -(v  +  a)  (4) 
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Figure  2.  Nomenclature  Used  in  the  Taylor  Analysis 

The  term,  1,  describes  the  rate  of  foreshortening  of  the 
undeformed  section  of  the  rod  and  can  be  written 

i  =  "(V  +  X)  (5) 

By  applying  Newton's  second  law  to  the  undeformed  portion  of  the 
rod  the  equation  of  motion  can  be  written  as 

dv  =  -  Y  (6) 

dt  (pi) 

Where  Y  and  p  are,  respectively,  the  material  yield  stress 
and  mass  density. 

Taylor  continues  the  analysis  by  writing  equations  describing 
conservation  of  mass  and  momentum  across  the  plastic  wave  front. 

A  differential  slice  of  the  undeformed  portion  of  the  rod,  dX, 
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Figure  3,  with  cross  sectional  area  A0,  crosses  the  plastic  front 
and  is  now  contained  in  the  volume  described  by  the  new  area,  A, 
and  the  differential  thickness,  dh.  The  elemental  length,  dX, 
can  be  written  in  terms  of  the  undeformed  section  using  the 
relationship,  dX  =  -di ,  Figure  2(a).  The  equation  for  the 
conservation  of  mass  can  be  written  as 

pAdh  =  pA0dX  (7) 

Dividing  both  sides  by  pdt,  Equation  (7)  becomes 

Adh  =  A0dX  (8) 

dt  dt 

After  substituting  for  dX  in  terms  of  dl ,  Equation  (8)  can  be 
rewritten  as 

AX  =  -  A0(i)  (9) 

Substituting  from  Equation  (5)  gives 

AX  =  A0(v  +  X)  (10) 

Taylor  assumes  the  material  behaves  in  such  a  way  that  when 
it  crosses  the  plastic-  front  it  comes  to  rest  instantaneously.  • 
This  assumption  imposes  a  condition  on  the  model  that  describes 
the  intermediate  states  of  the  event  as  having  a  strain 
discontinuity  at  the  elastic/plastic  interface,  Figures  1 
through  3.  This  strain  discontinuity  is  created  by  the 


Figure  3.  Schematic  Diagram  Showing  the  Volume  of  Material  That 
Passes  into  the  Plastic  Zone 
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instantaneous  change  in  cross  sectional  area  of  the  material 
passing  through  the  plastic  front. 

The  linear  momentum  equation  is  written 

pAvdSL  =  Y  ( A  -  A0)dt  (11) 

The  left  hand  side  of  the  equation  is  the  change  of  momentum  of 
the  differential  element,  dil,  having  an  initial  velocity,  v,  and 
a  final  velocity  of  zero.  The  right  hand  side  is  the  impulse 
term,  where  the  force  is  calculated  from  the  stress  in  the  body, 
assumed  uniform  over  the  cross  section,  times  the  relative 
change  in  area. 

To  construct  an  analysis  based  on  the  postmortem 
measurements  of  the  yield  boundary,  Taylor  must  describe  the 
motion  of  the  elastic/plastic  interface.  By  combining  Equations 
(6),  (10),  and  (11),  Taylor  showed  that  the  plastic  front  moves 

approximately  linearly  with  time.  Having  determined  this  from 
the  analysis,  he  subsequently  imposes  this  as  a  constraint  on 
the  model  when  developing  the  expression  for  the  yield  stress. 

The  expression  Taylor  uses  for  computing  the  yield  stress  is 
generated  from  the  equation  of  motion  of  the  undeformed  portion 
of  the  rod,  Equation  (6) .  The  independent  variable,  however, 
has  been  changed  from  that  of  time  to  the  incremental  change  in 
the  length  of  the  undeformed  section  of  the  rod.  This  gives 

dv  =  dv  dj.  =-  Y  (12) 

dt  dfl.  dt  (/si) 

Equation  (5)  can  be  substituted  into  Equation  (12)  to  give 

dv  =  Y  (13) 

dfl.  pfl (v  +  \ ) 

After  separating  the  variables,  integrating  and  substituting  the 
appropriate  initial  and  final  conditions,  the  equation  becomes 

Y_  inf  Lf  1  =  -  JL_  V2  -  VX  (14) 

p  L — E — J  2 
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where  Lf  is  the  final  length  of  the  specimen. 


To  eliminate  the  constant  plastic  wave  speed,  \,  Taylor 
assumes  that  the  rear  end  deceleration  of  the  rod  is  also 
constant.  This  assumption  allows  the  duration  of  the  ialpuCt 
event,  T,  to  be  calculated  two  ways.  First,  from  the  assumption 
of  constant  plastic  wave  speed 

T  =  (Lf  -  Af)  (15) 

_ 

Here  the  term,  If,  is  the  final  length  of  the  undeformed  segment 
of  the  specimen.  Then  from  the  assumption  of  constant 
deceleration 

T  »  2(L  -  Lf)  (16) 

V 

Eliminating  T  between  Equations  (15)  and  (16)  the  plastic  wave 
speed  is  determined  in  terms  of  the  impact  velocity  and  the  final 
specimen  geometry 

X_  «  (Lf  -  lt)  (17) 

V  2jZ  -  if) 

Using  Equation  (17)  to  eliminate  the  plastic  wave  speed  in 
Equation  (14),  the  yield  stress  is  determined  as  a  function  of 
density,  impact  velocity,  original  length,  final  undeformed 
length,  and  the  final  total  length.  This  gives  the  following 
expression  for  the  yield  stress: 

Y  =  pV*  (L  -  Lf)  1  (18) 

2  (L  -  Lf)  lnr  r  "I 

Nfj 

3.  IMPROVED  ANALYSIS  BY  TAYLOR 

In  a  second  approach,  Taylor  concedes  that  the  rear  end 
deceleration  is  not  constant.  He  suggests  a  more  exact  measure 
of  the  flow  stress  can  be  made  by  applying  a  correction  factor 
to  the  values  determined  in  the  above  simplified  analysis.  This 
correction  factor  is  calculated  from  the  error  introduced  by 
assuming  that  the  rear  end  deceleration  was  uniform. 
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To  determine  the  correction  factor,  Taylor  establishes  a  set 
of  equations  relating  the  plastic  wave  speed  and  the  rear  end 
motion  to  the  length  of  time  for  the  deformation  to  occur.  This 
set  of  expressions  can  be  solved  when  1 f/L  and  Lf/L  are  known 
quantities.  To  simplify  the  process,  a  graph  was  constructed 
with  Hf/L  as  the  ordinate  and  Lf/L  as  the  abscissa,  Figure  4, 
Reference  1.  The  appropriate  correction  factor,  Y/Yf,  can  be 
quickly  determined  from  the  coordinate  position  on  the  graph. 

The  term  Yf  describes  the  yield  stress  calculated  by  Equation 
(18).  In  general,  the  more  exact  analysis  will  increase  the 
yield  stress  level  of  the  material. 

Taylor's  analysis  predicts  that  the  cross  sectional  area  of 
the  deformed  region  will  vary  in  a  uniform  manner.  At  the 
elastic/plastic  interface,  where  the  state  of  stress  exceeds  the 
yield  strength  of  the  material,  the  material  will  deform 
radially  an  amount  which  is  dependent  on  the  current  velocity  of 
the  undeformed  segment  of  the  specimen.  Taylor  imposes  the 
condition  that  the  deceleration  of  the-  rear  end  be  uniform; 
therefore,  the  cross  sectional  area  in  the  deformed  region  must 
be  changing  uniformly.  In  reality,  the  specimen  profile  in  some 
materials  will  be  noticeably  nonuniform  and  will  depend  greatly 
on  the  strength  of  the  material  and  the  velocity  of  the  impact. 
Such  discrepancies  between  predicted  deformation  geometry  and 
actual  observation  is  related  to  the  assumptions  regarding  the 
material's  stress-strain  relationship  and  to  the  effect  of 
radial  inertia.  A  rigid,  plastic  material  behavior  model 
neglects  the  effects  of  complex  material  behavior  at  high 
strain  rates.  A  more  comprehensive  constitutive  model  might 
contain  terms  to  predict  such  phenomena  as  strain  hardening, 
thermal  softening,  strain  rate  effects,  dislocation  motion, 
grain  size  effects,  and  possibly  the  effects  due  to  the 
formation  of  deformation  twins. 
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4.  HAWKYARD 'S  ANALYSIS 

Hawkyard,  Reference  7,  modified  Taylor's  analysis  of  the 
problem  by  employing  an  energy  balance  expression  for  a 
differential  element  crossing  into  the  plastic  zone.  The  work 
comes  from  slowing  the  elemental  slice  of  the  specimen  from  its 
current  velocity  v  to  zero.  Additional  work  is  supplied  by  the 
force  acting  on  the  elemental  piece  by  the  undeformed  section. 
These  two  terms  combine  to  form  the  following  expression  for  the 
energy : 

Enerav  -  -Ac 1 o v2  -  YA0(-fl.-X)  (19) 

2 


Equation  (4)  can  be  used  to  make  substitutions  in  both  terms  to 
give 

Energy  =  Ac  ( v  \  )  ov2  +  YA0(v)  (20) 

2 


Hawkyard  equates  this  expression  to  the  rate 
work  in  uniaxial  compression  on  an  element. 


An ( v  -  X ) Yin 


oJ 


A^fv  *  X  )  p  v* 


of  doing 
-  YA0(v) 


plastic 


(21) 


Using  mass  conservation  across  the  plastic  front  and  force 
equilibrium  of  the  rigid  segment,  an  expression  can  be  developed 
that  predicts  the  profile  geometry.  Hawkyard' s  analysis  has 
achieved  a  certain  measure  of  success  in  predicting  a  concave 
profile,  which  has  been  observed  in  experiments,  Figure  5, 
Reference  7. 
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Figure  5.  Profile  Geometry  Predicted  By 

Hawkyard ' s  Energy  Balance  Analysis 


5.  a/ 1 B  ANALYSIS 

Critiques  of  the  analysis  by  Taylor  suggest  that  the 
momentum  equation  across  the  elastic  plastic  interface  was  not 
an  accurate  application  of  the  conservation  of  momentum.  Jones, 
Gillis,  and  Foster,  Reference  8,  report  that  a  more  precise 
expression  for  the  linear  impulse/momentum  equation  should 
contain  a  term  accounting  for  the  mass  loss  in  the  undeformed 
portion  of  the  rod,  Figure  6,  Reference  8.  At  some  time,  t, 
after  impact  the  undeformed  length  can  be  given  by 

A  =  L  -  X  (22) 

At  a  small  increment  of  time  later,  t+At,  the  undeformed  segment 
has  shortened  by  an  amount  equal  to  AX.  If  the  physical  system 
is  defined  to  be  the  initially  undeformed  portion  of  the  rod, 
the  impulse  equation  gives 

pA0AXu  +  pA0(L  -  X  -  AX) (v  +  AV)  -  pA0(L  -  X)v 

=  (  2  P  ±  A  P )  A  t 
2 
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(23) 


The  first  term  cn  the  left  side  is  the  momentum  of  the 
differential  element,  iX,  at  a  velocity  u,  as  it  passes  into  the 
plastic  zone.  The  second  term  accounts  for  the  momentum  of  the 
undeformed  segment,  which  has  undergone  a  mass  transfer,  and  has 
a  new  velocity  v  +  iv.  The  mass,  transferred  from  the 
undeformed  segment,  is  defined  in  terms  of  the  original  cross 
sectional  area,  A0,  and  the  differential  change  in  the  extent  of 
the  plastic  zone,  aX,  relative  to  the  original  configuration. 

v 


v  + A  v 


p  +  AP 

• 

— 

- -  L-X  -  AX - - 

Figure  6.  Mass  Transfer  From  the  Undeformed 
Section  of  the  Specimen 

The  last  term  on  the  left  is  the  original  momentum  of  the  system. 
The  right  hand  side  of  the  equation  is  the  impulse  applied  during 
the  time  increment. 

A  new  form  of  the  equation  of  motion  of  the  undeformed 
segment  of  the  rod  results  from  Equation  (23).  After  dividing 
both  sides  by  At,  taking  limits,  and  making  appropriate 
substitutions  Equation  (23)  becomes 

Hv  +  i (v  —  u)  =  P  (24) 

P  Aq 


i 
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The  force,  P,  can  be  replaced  by 


P  =  a  A  =  jA«  (25) 

(1  - ~ e ) 

Here  a  is  the  engineering  stress;  and  e  is  the 
longitudinal  compressive  strain.  Now 

1 V  +  0.  ( v  -  U)  =  a  ( 2  6) 

p ( 1  +  e) 

describes  the  motion  of  the  rigid  portion  of  the  specimen. 

If  there  is  no  change  in  particle  velocity  across  che 
interface,  v=u,  and  the  longitudinal  compressive  strain,  e,  is 
set  equal  to  zero,  then  Equation  (26)  reduces  to  Equation  (6) 
used  in  Taylor's  analysis.  From  the  conservation  of  mass, 
Equation  (10),  it  is  evident  that  Taylor  intended  the  model  to 
describe  the  elastic/plastic  front  as  a  strain  discontinuity. 

This  led  Jones  et  al,  Reference  8,  to  criticize  the  methods  used 
by  previous  investigators.  It  is  the  velocity  discontinuity,  v  * 
u,  across  the  interface  that  causes  the  deformation. 

Jones  et  al  developed  the  analysis  by  assuming,  as  Taylor 
did,  that  the  material  crossing  from  the  rigid  portion  of  the 
rod  into  the  plastic  zone  is  brought  instantaneously  to  rest, 
u=0 .  The  equation  of  motion  can  then  be  written  as 

dflv)  =  <7  (27) 

dt  p(l  +  e) 

Using  a  Lagrangian  coordinate  system,  the  strain  at  the  plastic 
front  can  be  described  by  a  relationship  among  the  kinematic 
variables . 

el  =  (v  -  u)  (28) 


Using  the  assumption  that  the  material  in  the  plastic  zone  has 
zero  velocity,  u=0 ,  the  strain  can  be  described  in  terms  of  the 
rear  end  velocity  divided  by  the  rate  of  foreshortening  of  the 
rigid  portion  of  the  rod.  If  the  plastic  wave  speed  is  assumed 
constant,  the  rate  of  foreshortening  can  be  described  by 
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Equation  (5).  The  strain  can  now  be  described  by  the  rear  end 
velocity  and  the  plastic  wave  speed.  Substituting  for  the  strain 
in  the  equation  of  motion  gives 


dULv).  =  <7  (v  ±.A1  (29) 

dt  p\ 

Using  the  chain  rule,  Equation  (29)  can  be  written  making  dl  the 
independent  variable 

d(lv)  =  dflv)  dl  =  a  (v  +  X )  (30) 

dt  dl  dt  pX 

Which  can  be  further  manipulated  after  substituting  from 
Equation  (5)  to  give 

dllvj  =  Z2.  (31) 

dl  p  x 

Following  the  use  of  the  product  rule,  Equation  (31)  can  be 
written  as 

Idv  +  v  =  -a  (32) 

dl  pX 

The  variables  can  now  be  separated,  and  the  resulting 
expression  integrated,  to  give  an  explicit  expression  for  1  in 
terms  of  the  velocity,  v. 


V  pX 


Before  the  right  hand  side  of  Equation  (33)  can  be 
integrated  a  constitutive  relationship  must  be  invoked  to 
specify  a.  Using  the  perfectly  plastic  material  model  gives 

a  =  -  Y  (34) 

Jones  et  al,  make  the  assumption  that  the  dynamic  yield  stress 
is  approximately  a  constant  for  an>  given  impact  velocity. 

After  substituting  Equation  (34)  into  Equation  (33),  the 
undeformed  length  is  given  as  a  function  of  the  impact  velocity, 
the  density,  the  yield  strength,  and  the  plastic  wave  speed, 
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V 


where 

c2  =  _Y 

p 

Evaluating  Equation  (35)  gives 


L  =  X  -  V  (36) 

J i  c2  -  v 

X 


After  substituting  the  appropriate  final  conditions,  Equation 
(36)  can  be  manipulated  into  the  form 


where 


(37) 


The  parameters,  a  and  £,  are  the  nondimensional  forms  of  the 
yield  strength  and  the  plastic  wave  speed,  respectively. 
Equation  (37)  contains  the  unknowns,  X  and  Y,  and,  therefore,  a 
second  independent  equation  is  required. 


The  second  equation  can  be  generated  from  the  equation  of 
motion  used  in  the  form  given  by  Equation  (29) .  After  using  the 
product  rule,  Equation  (29)  can  be  rewritten  as 

Idv  +  vdfl.  =  o  ( v  ±  X )  (38) 

dt  dt  p X 

Substituting  from  Equations  (5)  and  (34)  into  Equation  (38),  then 
separating  the  variables  gives 
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V 


tf  - 


_ Hdv 

Y (v  +  X )  +  v(v 

pK 

0 


x) 


(39) 


Where  the  limits  have  been  chosen  to  reflect  the  final 
conditions,  v=0,  and  t=tf.  The  velocity  dependent  term,  1,  can 
be  expressed  by  Equation  (36)  and  substituted  into  Equation  (39) 
to  give 


V 

r  dv 

Lf  Y  -VI 

(v  +  x) 

'  Y  -r  V' 

J  j 

^pX 

0 


After  evaluating  the  integral,  making  appropriate  substitutions, 
and  performing  some  algebraic  manipulation,  Equation  (40)  gives 


Lf  -  If 


L 


0(1  -  &f/L)2  ' 

.(1  +  0(1  -  lf/L). 


0(Af/L)  (1  -  Af/L) 

In 

0  ( A  f/L) 

[  (i  +  0(1  -  af/L))j 

l  1  +  0  J 

Using  suitable  experimental  data,  Equation  (41)  can  be  solved 
numerically  for  the  parameter  0.  By  substituting  the  value  of  0 
into  Equation  (37)  the  parameter  a  is  determined.  Knowing  a  and 
0,  the  material  flow  stress,  Y,  and  the  Eulerian  plastic  wave 
speed,  X,  are  easily  found. 


6 .  EXPERIMENTAL  MATRIX 

An  experimental  matrix  was  constructed  to  test  the  model 
proposed  by  investigators  Jones,  Gillis,  and  Foster,  Reference  8. 
This  matrix  consisted  of  several  types  of  materials.  The 
specimen  geometries  used  ih  this  investigation  covered  a  wide 
range  of  aspect  ratios  (length/diameter) .  The  influence  of 
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specimen  geometry  is  a  factor  that  has  been  debated  in  the  liter¬ 
ature  by  several  investigators,  References  9  and  10. 

The  materials  chosen  consisted  of  two  types  of  aluminum  and 
copper.  The  aluminum  materials  were  alloys  2024-T4  and  6061-T6. 
The  copper  materials  were  Oxygen  Free  High  Conductivity  (OFHC) 
copper  and  a  Phosphorous  Deoxidized  Tellurium  copper  alloy 
( DPTE ) ;  both  coppers  were  received  in  a  half -hard  condition.  The 
materials  were  chosen  with  regard  to  the  sensitivity  of  the 
yield  stress  to  the  rate  of  straining.  The  6061-T6  material  is 
known  to  be  relatively  insensitive  to  rate  effects,  whereas  the 
2024-T4  and  the  OFHC  copper  are  known  to  be  rate-sensitive 
materials.  The  DPTE  copper  is  a  solution  strengthened 
copper/tellurium  alloy,  chosen  for  comparative  purposes  with  the 
higher  purity  copper. 

Mechanical  properties  of  the  test  materials  were  determined 
from  tensile  tests  and  hardness  measurements.  A  typical  load 
versus  time  chart  produced  from  a  tensile  test  of  each  material 
can  be  found  in  Appendix  D.  Results  of  these  test  and  hardness 
measurements  are  reported  in  Table  1. 

TABLE  1.  MECHANICAL  PROPERTIES  (QUASI-STATIC) 


MATERIAL 

YIELD 

STRENGTH 

MPa 

TENSILE 

STRENGTH 

MPa 

PERCENT 

ELONGATION 

HARDNESS 

Rb 

OFHC  COPPER 

— 

350 

17.6 

33.5 

DPTE  COOPER 

— 

300 

12.2 

32.0 

2024-T4 

400 

500 

23.0 

75.5 

6061-T6 

315 

340 

18.1 

45.7 
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SECTION  III 


EXPERIMENTAL  APPARATUS 


1 .  BACKGROUND 

In  order  to  conduct  the  proposed  test  matrix,  a  suitable 
experimental  apparatus  was  required.  Several  factors  determined 
by  experience  with  a  compressed  air  gun  made  it  imperative  that 
a  new  apparatus  be  constructed.  The  baseline  operational 
requirements  called  for  a  flexible,  repeatable,  and  efficient 
system  design.  The  perspective  from  which  the  new  apparatus  was 
developed  can  be  better  understood  following  a  brief  summary  of 
the  configuration  and  capabilities  of  the  compressed  air  gun 
system. 

The  gas  gun  system  was  simply  a  high  pressure  tank,  or 
vessel,  attached  to  the  breech  end  of  a  gun  barrel.  The  muzzle 
end  of  the  gun  barrel  was  permanently  fixed  against  the  side  of 
a  holding  tank.  The  target  was  positioned  in  the  holding  tank 
at  a  standoff  distance  of  approximately  38  cm  from  the  muzzle. 

The  operation  of  the  gas  gun  required  several  steps.  First, 
a  thin  metal  diaphragm  (bursting  disk)  was  placed  over  the 
discharge  orifice  of  the  pressure  tank.  Next,  a  specimen  was 
loaded  into  the  gun  barrel,  after  which  the  barrel  and  the 
pressure  vessel  were  bolted  together,  sealing  the  pressure  tank 
orifice.  At  this  point,  the  vessel  could  be  pressurized  to  the 
desired  level.  Once  the  prescribed  pressure  level  was 
established,  the  diaphragm  was  punctured,  via  a  mechanical 
striker,  allowing  the  pressurized  gas  to  escape  and  accelerate 
the  specimen  through  the  barrel. 

A  number  of  deficiencies  were  experienced  with  the  use  of 
such  an  apparatus.  At  best,  the  gas  gun  system  efficiency  was 
quite  poor.  After  each  shot,  it  was  necessary  to  unbolt  the 
pressure  tank  from  the  barrel.  In  addition,  it  took  several 
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minutes  to  build  up  the  operating  pressure  needed  for  the  test. 
In  terms  of  operational  capability,  the  system  had  very  poor 
repeatability . 

The  lack  of  repeatability  was  caused  by  two  main  problems: 
the  bursting  disk  often  ruptured  prematurely  and  variations  in 
the  machining  of  the  specimen  allowed  gas  to  escape  past  the 
specimen.  For  shots  made  at  equal  tank  pressure,  the  impact 
velocities  would  vary  widely  if  leakage  occurred. 

The  apparatus  had  built-in  limitations  with  regard  to 
flexibility.  The  pressure  vessel  had,  for  safety  reasons,  an 
upper  limit  of  5.17  MPa.  Hence,  the  total  available  energy  was 
fixed.  For  high  density  materials,  this  energy  level  was  not 
sufficient  for  experiments  in  the  desired  velocity  range. 

Often,  the  lack  of  adequate  energy  was  manifested  by  oblique 
impact  with  the  target.  Since  the  muzzle  and  target  positions 
were  fixed,  the  standoff  distance  could  not  be  adjusted.  For 
these  reasons,  a  new  apparatus,  Figure  7,  was  developed  based  on 
the  use  of  a  smokeless  gun  powder  as  a  propellant. 

2.  DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS 

The  launch  tube  was  machined  from  4340  steel  and  was  center 
bored  to  an  inside  diameter  of  7.620  mm.  Two  ports  were  drilled 
into  the  gun  barrel  near  the  muzzle.  The  ports  were  2.54  cm 
apart  and  were  tapped  to  receive  pressure  transducers. 

Monitoring  the  electrical  signals  of  the  transducers  was  one 
method  of  determining  the  projectile  velocity  near  the  muzzle  of 
the  launch  tube. 

The  breech  end  of  the  launch  tube  was  chambered  to  accept  a 
0.308  caliber  cartridge  case.  The  bullet  and  original  powder 
charge  were  removed  from  a  0.308  rifle  cartridge  so  that  the 
primed  cartridge  case  could  be  used.  In  place  of  the  original 
powder  charge,  a  specific  quantity  of  smokeless  powder  (Red  Dot) 
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Figure  7 .  A  Taylor  Anvil  Test  Apparatus 

• 

was  loaded.  The  propellant  was  covered  with  !a  small  wad  of  cot¬ 
ton.  The  cotton  was  packed  against  the  powder  and  primer  of  the 
cartridge.  The  function  pf  the  cotton  was  to  ensure  a  uniform 
burn  rate  of  the  propellant  from  shot  to  shot  by  keeping  the 
powder  charge  in  place  against  the  primer.  The  powder  charges 
used  for  this  work  varied  from  1  to  5  grains  (1  gram  =  15.4 
grains).  By  comparison,  the  normal  0.308  rifle  cartridge  powder 
charge  is  approximately  40  grains.  Consequently,  a  large  portion 
of  the  volume  was  tilled  by  the  cotton. 

To  lock  the  cartridge  in  the  chamber,  the  launch  tube  had 
been  externally  threaded  to  facilitate  the  mounting  of  an  end  cap 
on  the  breech.  The  end  cap  contained  a  through  hole  for  the 
firing  pin  in  order  to  make  contact  with  the  cartridge,  Figure  8. 
Actuation  of  the  firing  pin  was  by  means  of  an  electric  solenoid, 
coi. trolled  from  the  firing  room. 
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To  prevent  the  expanding  powder  gas  from  leaking  past  the 
specimen,  a  plastic  obturator  was  used  to  form  a  seal.  The  obtu¬ 
rator  was  positioned  between  the  cartridge  and  the  specimen, 
Figure  8.  The  end  of  the  obturator  nearest  the  cartridge  was 
hollowed  to  facilitate  radial  expansion  of  the  remaining  material 
under  pressure.  Thus,  the  expanding  gas  deforms  the  obturator  to 
form  a  seal  against  the  gun  bore  wall.  In  this  manner,  the 
energy  of  the  propellant  was  used  entirely  for  accelerating  the 
specimen. 


Figure  8.  Firing  Line  Components 

The  launch  tube  was  fixed  in  position  by  a  set  of  v-block 
mounts.  This  construction  allowed  rapid  change  of  the  standoff 
distance  between  the  muzzle  and  the  target.  Typically,  the 
standoff  was  in  the  range  of  7  to  20  cm. 

The  target  was  a  23-cm  diameter,  by  20  cm  in  length,  cylinder 
of  hardened  4340  steel.  Both  ends  were  machined  parallel  and  lap 
finished.  The  design  of  the  target  and  target  rest,  Figure  9, 


optimized  the  available  surface  area  for  impact  tests.  After 
each  test,  the  target  could  be  rotated  to  provide  a  new  surface 
for  impact.  If  no  permanent  deformation  occurred  in  the  target 
surface,  several  tests  were  conducted  using  the  same  area  of  the 
anvil.  After  completing  one  revolution  of  the  target,  the 
center  line  of  the  cylindrical  anvil  could  be  lowered  with 
respect  to  the  projectile  flight  line  to  provide  a  new  surface 
area  for  impact.  To  lower  the  target,  a  layer  of  the  polypropo- 
lux  base  material  was  removed  from  under  the  target  rest.  This 
process  of  lowering  the  anvil  center  line  could  be  continued 
until  the  entire  surface  area  was  used.  Subsequently,  the  paral¬ 
lel  faces  of  the  anvil  could  be  reversed  and  the  process  repeated 
prior  to  remachining.  This  optimized  the  available  area  on  the 
target  and  reduced  the  down  time  required  for  remachining  and 
polishing  of  the  surface. 


Figure  9.  Target  Design  and  Photographic  Configuration 

The  target  and  the  muzzle  of  the  launcher  were  covered  by  an 
aluminum  housing.  Two  slots  were  cut  into  the  housing  parallel 
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to  the  flight  line  of  the  projectile.  These  slots  were  covered 
with  6.3  5  mm  thick  plexiglass,  providing  w’indows  vhiwh  allowed 
the  incoming  projectile  and  its  subsequent  deformation  to  be 
photographed.  The  housing  prevented  the  rebounding  projectile 
from  causing  undesirable  damage.  The  inside  of  the  housing  was 
lined  to  prevent  secondary  deformation  from  occurring  on  the 
specimen.  However,  a  large  number  of  specimens  were  slightly 
deformed  on  the  rear  end  of  the  projectile  because  of  impact 
with  the  muzzle  face  of  the  barrel  after  rebounding  from  the 
target. 

3.  VELOCITY  MEASUREMENT  TECHNIQUES 

During  the  test,  as  many  as  three  techniques  were  used  to 
determine  the  projectile  velocity.  As  previously  mentioned,  the 
gun  barrel  was  instrumented  with  two  piezoelectric  pressure 
transducers.  The  outputs  from  the  transducer  amplifiers  were 
fed  into  a  dual  trace  oscilloscope.  From  the  oscilloscope,  the 
time  required  for  the  expanding  gas  of  the  propellant  to  pass 
between  the  two  transducers  could  be  measured.  Knowing  the 
time,  and  the  distance  between  the  transducer  ports,  the 
velocity  was  easily  calculated. 

A  learning  period  was  necessary  to  determine  what  threshold 
sensitivity  level  to  set  on  the  transducer  amplifiers  for  the 
various  pressure  levels  seen  in  the  barrel.  If  the  amplifiers 
were  set  too  sensitively,  the  signal  would  be  erratic,  possibly 
from  the  elastic  wave  in  the  gun  itself.  If  the  sensitivity  was 
set  too  low,  no  response  would  be  obtained.  In  both  cases,  the 
necessary  sensitivity  level  was  always  pressure  dependent. 
Eventually,  the  amount  of  propellant  became  the  best  source  of  a 
priori  information  on  how  to  set  the  sensitivity  level. 

Error  in  the  velocity  found  by  using  the  pressure 
transducers  can  be  analyzed  by  examining  the  source  for  the 
data.  The  signals  from  the  pressure  transducers,  located  near 
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the  muzzle  end  of  the  launch  tube,  are  sent  through  a  signal 
amplifier  before  being  displayed  on  a  digital  oscilloscope, 

Figure  10.  The  signals  displayed  on  the  oscilloscope  can  be 
measured  with  the  sweep  cursor  to  determine  a  time  interval 
necessary  for  calculating  the  muzzle  velocity,  Figure  11(a).  It 
was  considered  for  this  analysis  that  uncertainties,  or 
systematic  errors,  in  the  measuremert  system,  i.e.,  line  noise, 
transducer  and  electronic  circuit  response  characteristics,  etc., 
were  small  by  comparison  to  the  error  introduced  by  misidentify- 
ing  the  starting  and  ending  points  of  the  interval  to  be  meas¬ 
ured. 

The  Nicolet  oscilloscope  displays  a  digitized  form  of  the 
s  .gnal  as  a  series  of  individual  points.  Each  point  represents 
a  segment  of  time,  which  for  this  study  was  1  fis.  The  total 
time  displayed  could  vary  depending  on  the  necessary  resolution. 

A  desirable  wave  form  would  show  a  smooth  baseline,  before  the 
obturator  passes  the  transducer  port,  followed  by  a  sharp  point, 
the  response  of  the  transducer  to  the  propellant  gas  pressure, 
Figure  11(a).  A  poor  signal  would  have  a  slowly  changing 
response,  or  curved  wave  form,  prior  to  a  rapid  deflection, 

Figure  11(b). 

The  former  type  of  signal  would  provide  the  necessary  infor¬ 
mation  to  determine  the  time  ir.terv*- 1  _o  2  ms.  The  accu¬ 

racy  from  a  measurement  of  a  projectile  having  a  velocity  of  200 
m/s,  using  this  type  of  wave  form  would  be  ±3  m/s.  The  second 
type  of  wave  form  would  have  greater  uncertainty  because  of  the 
subjectivity  of  establishing  a  suitable  baseline  and  break  point 
for  the  signal. 

The  second  signal  from  the  amplifier  was  sent  to  a  multitrack 
recorder.  The  transducer  signals  were  played  back  to  provide  an 
additional  method  of  digitizing  the  data.  A  strip  chart  output 
of  the  recorded  signals  and  a  reference  timing  signal  was  pro¬ 
duced.  These  data  were  digitized  by  establishing  a  timing  base- 
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Figure  10.  Block  Diagram  of  the  Signal  Flow  for  the 
Pressure  Transducer  Measurement  System 


#1  PRESSURE  TRANSDUCER  *1  PRESSURE  TRANSDUCER 


(a) 


(b) 


Figure  11.  Typical  Wave  Forms  Generated 
by  the  Pressure  Transducers 
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line  from  the  reference  signal  (100  kHz  -1Hz).  The  starting  and 
ending  points  were  identified  and  the  interval  of  time  could  be 
established.  Similar  to  the  digital  oscilloscope,  the  choice  of 
starting  and  ending  points  on  the  interval  was  subjective.  The 
quantitative  amount  of  error  possible  was  dependent  on  the  signal 
wave  form.  In  general,  measurements  taken  directly  from  the 
digital  oscilloscope  were  considered  to  be  the  more  accurate  of 
the  two  methods,  primarily  because  of  the  higher  level  of  resolu¬ 
tion  of  the  wave  form. 

As  a  second  measure  of  the  projectile  velocity,  a  high  speed 
movie  camera  was  also  used.  The  16  mm  camera  was  capable  of 
10,000  frames  per  second.  Using  a  1/4  frame  format,  a  frame 
rate  of  40,000/s  could  be  obtained.  A  shadowgraphic  technique 
was  used  to  photograph  the  incoming  projectile  and  its 
deformation  at  the  anvil  face.  To  do  so  required  a  light  source 
located  behind  the  specimen,  as  shown  in  Figure  9. 

A  fiducial  marker  (cylindrical  magnet) ,  of  known  diameter 
and  length,  was  placed  on  the  anvil  face  directly  above  the 
point  of  impact.  With  this  in  the  field  of  view  of  the  camera, 
it  was  possible  to  determine  the  degree  of  obliquity  the  camera 
had  with  the  anvil  face.  By  through-the-lens  alignment,  the 
camera  was  set  as  near  parallel  as  possible  to  the  anvil  face. 

As  a  check,  the  profile  image  of  the  fiducial  recorded  on  the 
film  could  be  measured  to  determine  if  the  length,  viewed  from 
the  position  of  the  camera,  has  been  either  elongated  or  shor¬ 
tened.  The  actual  position  of  the  anvil  face,  at  the  flight 
line,  could  then  be  located. 

To  determine  a  velocity  from  the  film,  it  was  necessary  to 
find  the  time  for  the  particular  distance  traversed.  The  time 
was  calculated  by  multiplying  the  number  of  frames  by  the  frame 
rate.  An  average  value  of  the  frame  rate  was  determined  by 
timing  marks  recorded  on  the  film.  The  distance  traveled  could 
be  measured  using  a  film  analyzer.  A  reference  point  at  a 
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particular  frame  was  established  by  positioning  the  vertical 
cross  hair  on  the  leading  edge  of  the  specimen,  then  setting  the 
digital  counter  to  zero.  By  moving  the  cross  hair  to  the  same 
location  on  the  specimen  at  a  different  frame,  a  distance 
traversed  could  be  measured.  A  scaling  factor,  for  the 
magnification,  was  determined  by  measuring  the  specimen  diameter 
in  digital  counter  units.  Multiplying  the  scaling  factor  with 
the  counter  distance  measured  for  the  traverse  distance  gives 
the  actual  distance. 

Limits  on  the  accuracy  of  such  a  technique  were  from  two 
sources.  First,  the  frame  rate  had  to  be  averaged  over  a  signif¬ 
icantly  larger  portion  of  film  than  was  used  in  the  velocity 
measurement.  Secondly,  the  location  of  the  cross  hair  on  the 
leading  edge  of  the  specimen  was  subjective.  The  motion  of  the 
specimen  during  the  actual  exposure  caused  the  image  to  be 
slightly  blurred.  However,  the  overall  technique  was  estimated 
to  be  accurate  to  within  ±10  m/s  over  the  range  of  velocities 
used  in  this  study. 

On  some  of  the  test  shots,  the  movie  camera  was  replaced 
with  a  high  speed  framing  camera  to  produce  a  detailed 
shadowgraphic  recording  of  the  impact  event.  This  high  speed 
Cordon  camera  was  operated  at  0.30  million  frames  per  second. 

At  this  framing  rate  the  resolution  of  the  data  generated  was 
one  photograph  every  3.3  fis.  By  comparison,  the  movie  camera 
produced  one  photograph  every  25  n s. 

Unfortunately,  the  complex  design  of  the  framing  camera 
allowed  only  82  frames,  35  mm  size,  to  be  recorded.  This  created 
a  timing  window,  160  ns,  in  which  several  events  had  to  occur. 

In  this  particular  camera  design,  when  it  came  up  to  the  desired 
operating  speed  the  shutter  was  automatically  triggered.  This 
created  the  situation  where  the  light  source  controlled  the  film 
exposure  level.  By  its  very  nature,  the  lighting  found  in  high 
speed  photography  requires  extremely  specialized  equipment.  For 
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this  camera,  it  was  necessary  to  produce  a  large  quantity  of 
light  energy  for  a  very  precise  period  of  time.  Too  little 
light  would  produce  poor  film  quality,  while  a  lengthy  exposure 
caused  overwriting  or  a  double  image. 

The  timing  window  also  required  that  the  light  source  be 
synchronized  with  the  projectile's  arrival  at  the  target  face. 
This  synchronization  was  accomplished  by  using  the  pressure 
transducers,  located  at  the  muzzle  of  the  launch  tube,  with  an 
appropriate  delay  circuit. 

A  third  measurement  system  for  the  velocity  relied  on 
infrared  beams  and  detectors  between  the  anvil  and  the  muzzle  to 
record  the  position  of  the  incoming  projectile.  The  hardware 
used  for  this  system  can  be  seen  in  Figure  7,  near  the  anvil 
face.  The  concept  of  the  system  was  to  use  the  beam  emitters  to 
produce  a  high  voltage  state  in  the  detector  circuit.  Once  the 
projectile  passed  into  the  beam,  the  detector  would  drop  to  a 
low  voltage  state.  Monitoring  the  voltage  states  of  the  two 
detectors  with  an  oscilloscope  provided  a  time  increment  for  the 
projectile  to  pass  between  the  beams.  Knowing  the  distance 
between  the  beams,  the  velocity  was  determined. 

In  this  system  the  limitations  were  related  to  the  response 
time  characteristics  of  the  circuitry.  To  reduce  the  error,  the 
components  used  were  individually  compared  with  their 
counterpart  to  ensure  both  detector  systems  responded  at  equal 
rates.  In  general,  the  velocities  measured  by  the  infrared 
detectors  were  considered  the  most  accurate  of  the  three 
methods.  All  of  the  methods  usually  gave  values  within  5 
percent  of  each  other  for  the  range  of  velocities  used. 

4 .  TEST  PROCEDURES 

The  basic  procedures  for  conducting  the  experiment  can  be 
broken  into  three  groups:  operations  prior  to  firing  the  gun, 
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the  actual  test,  and  recovery  and  postmortem  measurements.  A 
step  by  step  listing  can  be  found  in  Appendix  A. 

The  initial  operation  consisted  of  weighing  and  measuring 
the  diameter  and  length  of  the  specimen.  Knowing  the  mass  of 
the  specimen  and  the  velocity  of  interest,  a  propellant  charge 
could  then  be  specified.  Initially,  the  process  of  using  the 
smokeless  powder  was  one  of  trial  and  error.  However,  aft  r 
acquiring  sufficient  data  on  mass,  velocity,  and  propellant 
weight,  a  graph  was  produced  to  provide  a  quick  source  for  this 
information,  Figure  12. 

Prior  to  arming  the  launch  tube,  a  number  of  instrumentation 
checks  and  cleaning  operations  were  performed.  As  an  example, 
debris  often  fell  on  the  infrared  detector  lens,  which  had  to  be 
removed  for  the  device  to  operate  properly.  Other  operations 
included  rotating  the  anvil  and  checking  the  alignment  between 
its  face  and  the  launch  tube  muzzle.  At  this  point,  it  was 
appropriate  to  set  the  standoff  distance  between  the  muzzle  and 
the  target.  In  parallel  with  these  operations,  the  alignment  and 
loading  of  the  high  speed  movie  camera  were  normally  conducted. 
After  alignment  and  checkout  of  the  instrumentation,  the  aluminum 
housing  was  placed  over  the  target  and  the  muzzle  of  the  launch 
tube . 

The  actual  test  was  conducted  by  first  loading  the  specimen 
and  obturator  into  the  launch  tube.  A  gauge  was  used  to  position 
the  specimen  and  obturator  at  the  same  location  in  the  tube  for 
each  test.  The  cartridge  was  then  placed  into  the  chamber  behind 
the  obturator,  Figure  8.  The  end  cap  was  screwed  onto  the  launch 
tube  until  marks,  scribed  on  the  gun  and  cap,  were  aligned.  The 
electric  solenoid  with  the  firing  pin  was  then  placed  on  the  end 
cap.  Electrical  cables,  from  the  safe/arm  control  box,  were  then 
connected  to  the  solenoid. 
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Taylor  Experiment 

Propellant  Weight  vs  Velocity 

Propellant  Type:  Red  Dot  Smokeless  Powder 
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gure  12.  Propellant  Weight  Versus  Velocity  Curve  for  Different 
Specimen  Masses 


After  arming  the  device,  the  photo  lamps  were  turned  on.  At 
3  seconds  prior  to  firing,  the  movie  camera  was  switched  on. 

This  allowed  time  for  the  camera  to  come  up  to  its  maximum  fram¬ 
ing  rate  before  recording  the  impact  event.  At  firing,  an  elec¬ 
trical  signal  was  sent  to  the  solenoid,  which  then  drove  the 
firing  pin  into  the  primer.  The  detonation  of  the  primer  caused 
the  propellant  to  react,  producing  the  volume  of  expanding  gas 
necessary  to  accelerate  the  specimen.  Following  impact,  the 
specimen  was  recovered  from  under  the  aluminum  housing. 

The  material  used  in  this  investigation  was  received  as  rods 
8  mm  in  diameter  by  3.66  m  in  length.  Each  of  the  rods  were  cut 
into  5  smaller  sections.  From  each  section,  two  tensile 
specimens  of  the  material  were  machined,  a  total  of  ten  from 
each  rod.  The  remainder  of  the  section  was  machined,  on  a 
precision  lathe,  to  7.595  mm  in  diameter,  before  being  cut  into 
specimens  of  the  desired  lengths.  Table  2.  The  range  of  aspect 
ratios  covered  in  this  investigation  was  from  1.5  to  10. 


TABLE  2.  SPECIMEN  DIMENSIONS 
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5 .  POSTMORTEM  MEASUREMENTS 

For  making  postmortem  measurements,  a  gauge  was  developed  to 
provide  rapid  and  accurate  measurements  of  the  undeformed  length, 
Figure  13.  The  gauge  was  machined  from  a  piece  of  round  stock 
2.5  cm  diameter  by  6.3  cm  in  length.  A  7.595  mm  through  hole  was 
machined  in  the  gauge.  Then  it  was  counterbored  from  one  end 
until  only  a  very  thin  portion  of  the  original  diameter  hole  was 
remaining . 

Placing  the  undeformed  end  of  the  specimen  into  the  through 
hole,  measuring  7.595  mm,  the  plastic  zone  of  the  impact 
specimen  would  extend  out  of  the  top  of  the  gauge.  Taking  the 
difference  between  the  final  specimen  length  Lf  and  the  length 
of  the  plastic  zone  protruding  from  the  gauge  determines  the 
final  undeformed  specimen  length,  If.  This  method  proved  quite 
successful  at  reducing  scatter  in  the  data  when  compared  with 
previously  used  methods. 


SCALE 


SECTION  "A— A” 

Figure  13.  Detail  Drawing  of  the  Measurement  Gauge 


SECTION  IV 


ANALYSIS  OF  TEST  RESULTS 


1 .  INTRODUCTION 

The  primary  goal  of  the  experimental  work  was  to  provide  a 
aaua  case  to  us*=  with  various  computer  models.  Information 
contained  in  the  data  base  was  used  to  calculate  the  yield 
stress  and  to  construct  graphs  of  the  yield  stress  as  a  function 
of  various  experimental  parameters.  In  addition,  a  study  was 
made  of  the  sensitivity  of  the  computed  yield  stress  to 
perturbations  in  the  experimental  data.  From  this  study  an 
estimate  of  the  quantitative  amount  of  error  possible  in  the 
calculations  of  the  yield  stress  can  be  made.  This  report  will 
give  evidence  supporting  the  use  of  the  a/p  model,  Reference  8, 
for  predicting  the  material  behavior  under  impact  conditions. 

The  test  matrix  was  designed  to  generate  data  covering  a 
wide  range  of  velocities,  material  types,  and  specimen  aspect 
ratios.  The  experimental  data,  along  with  results  calculated 
from  the  various  analytical  models,  were  organized  by  material 
type  and  can  be  found  in  Appendix  B. 

The  range  of  impact  velocities  used  in  this  investigation 
was  120  to  330  m/s.  For  a  particular  material,  the  upper 
boundary  for  the  impact  velocity  was  limited  by  radial  cracking. 
For  pure  copper,  the  maximum  impact  velocity  was  approximately 
200  m/s.  The  DPTE  copper  alloy  had  a  maximum  velocity  of  160 
m/s.  For  the  aluminum  alloys,  the  2024-T4  material  had  a 
maximum  velocity  of  290  m/s,  while  the  6061-T6  had  a  maximum 
velocity  of  330  m/s. 
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2.  POSTMORTEM  MEASUREMENT  RELIABILITY 

The  reliability  of  the  experimental  measurements  can  be 
judged  by  monitoring  the  trends  in  the  data  as  a  function  of  a 
process  variable.  Figure  14  shows  how  the  nondimensional  forms 
of  the  final  length,  Lf/L,  and  the  final  undeformed  length, 

vary  with  the  impact  velocity  for  ^he  lour  test  materials. 
These  data  show  that  the  final  length  and  undeformed  length 
decrease  approximately  linearly  with  increasing  velocity  of 
impact . 

After  curve  fitting  the  data  by  linear  regression,  a  statis¬ 
tical  analysis,  Table  3,  reveals  the  level  of  uncertainty  to  be 
found  in  the  postmortem  measurements.  For  the  final  length 
measurement,  the  largest  average  deviation  and  maximum  deviation 
were  1.41  and  4.35  percent,  respectively.  A  similar  analysis 
revealed  a  considerable  amount  of  scatter  was  present  in  the  data 
for  the  final  undeformed  length.  The  largest  average  and  maximum 
deviations  were  6.75  and  42.09  percent,  respectively. 

The  primary  source  of  scatter  in  the  measurement  of  the 
final  undeformed  length  is  from  non-symmetric  deformation. 

Oblique  impact  with  the  target  was  often  found  to  be  the  source 
of  such  deformation.  When  a  specimen  impacts  the  target 
obliquely,  the  plastic  zone  will  extend  further  down  one  side  of 
the  specimen  than  on  the  other.  The  measurement  technique  for 
determining  the  final  undeformed  length  used  a  gauge  developed 
with  the  premise  that  the  specimen  deformed  in  an  axisymmetric 
manner.  Specimens  that  were  visually  identified  as  having 
impacted  obliquely  were  not  measured. 
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TABLE  3.  LINEAR  REGRESSION  AND  STATISTICAL  DATA 


MATERIAL 

DATA 

SLOPE 

(l/cm/s) 

Y-INTERCEPT 

(cm/cm) 

AVERAGE  % 
DEVIATION 

MAXIMUM  % 
DEVIATION 

OFHC 

Lf/L  vs  V 

-  0.0020 

1.124 

1.41 

4.35 

OFHC 

if/ L  vs  V 

-  G.CGlb 

0.644 

6.75 

42.09 

DPTE 

Lf/L  vs  V 

-  0.0022 

1.145 

0.81 

1.71 

DPTE 

lf/L  vs  V 

-  0.0032 

0.8327 

3.32 

7.22 

2024-T4 

Lf/L  vs  V 

-  0.0008 

1.088 

0.80 

2.48 

2024-T4 

lf/L  vs  V 

-  0.0014 

0.712 

3.70 

13.95 

6061-T6 

Lf/L  vs  V 

-  0.0012 

1.140 

0.83 

2.79 

6061-T6 

lf/L  vs  V 

-  0.0013 

0.712 

4.02 

10.36 

3.  STATIC  AND  DYNAMIC  STRENGTHS 

The  strength  of  the  materials,  at  low  strain  rates,  was 
determined  by  quasi-static  tensile  tests.  Specimens  for  the 
test  were  machined  from  various  sections  of  each  of  the  rods. 

The  specifications  for  the  fabrication  of  the  tensile  specimens 
can  be  found  in  Appendix  C.  The  tensile  tests  were  conducted  on 
a  screw-driven  Instron  machine  at  a  crosshead  displacement  rate 
of  8.5  x  10-3  mm/s.  Typical  load  versus  time  plots  can  be  found 
in  Appendix  D. 

The  data  obtained  from  the  impact  tests  was  used  in  a 
computer  program  which  determines  the  dynamic  yield  stress  of 
the  material.  The  program  was  designed  to  solve  for  the 
nondimens ional  plastic  wave  speed,  0,  and  the  nondimensional 
yield  stress,  a.  The  parameter  0  is  determined  by  Equation 
(41),  page  19.  The  computer  program  uses  an  interval  halving 
technique  to  solve  Equation  (41)  for  the  parameter  0.  The  par¬ 
ameter  a  is  determined  by  Equation  (37),  page  18,  and  is  directly 
calculable  once  0  has  been  found.  Knowing  a  and  0,  the  yield 
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stress  and  the  plastic  wave  speed  can  be  easily  determined  from 
their  definitions  given  following  Equation  (37),  page  18. 

The  values  of  the  yield  stress  from  tensile  test  and  those 
computed  by  the  a/p  model  are  shown  in  Table  4.  In  general,  the 
dynamic  yield  stress  is  approximately  1.5  to  2  times  the  quasi¬ 
static  value  of  the  yield  stress.  This  increase  in  the  yield 
strength  in  rhe  material  corresponds  to  a  change  in  the  strain 
rate  of  7-8  orders  of  magnitude. 

In  Table  4,  the  yields  stress  values  for  the  two  copper 
materials  were  not  reported.  The  response  of  the  copper 

TABLE  4.  COMPARISON  OF  STATIC  AND  DYNAMIC  YIELD 
STRESS  VALUES 


QUASI-STATIC 

( TENSILE  TEST) 

DYNAMIC 

(IMPACT  TEST) 

MATERIAL 

YIELD 

STRESS 

ULTIMATE 

STRESS 

YIELD 

STRESS 

OFHC 

— 

350  Mpa 

550  Mpa 

DPTE 

— 

300  Mpa 

520  Mpa 

2024-T4 

400  Mpa 

500  Mpa 

750-900  Mpa 

6061 -T6 

315  Mpa 

340  Mpa 

550  Mpa  | 

materials  to  loading  was  affected  by  the  high  dislocation 
density  pre-existing  in  the  material.  This  material  condition 
created  a  load  versus  time  curve,  Appendix  D,  in  which  the 
stress  level  rose  linearly  to  a  maximum  value  and  then 
immediately  decreased.  The  high  dislocation  density  in  the 
material  effectively  eliminated  the  work  hardening  region  of  the 
stress  versus  strain  curve. 

4.  SENSITIVITY  ANALYSIS 

An  estimate  of  the  error  found  in  the  computed  value  of  the 
yield  stress  can  be  determined  from  a  sensitivity  analysis. 
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This  analysis  was  conducted  by  estimating  the  level  of  accuracy 
associated  with  each  experimental  measurement,  e.g.,  the  final 
undeformed  length.  The  value  of  the  yield  stress  was 
recalculated  after  the  data  were  perturbed  an  amount  equal  to 
the  uncertainty  in  each  measurement.  A  quantitative  estimate  of 
the  possible  error  is  found  by  comparing  the  baseline  value  of 
the  yield  stress  with  those  found  in  the  sensitivity  analysis. 

The  measurements  of  specimen  geometry,  the  final  length,  and 
the  final  undeformed  length  were  estimated  to  be  accurate  to 
within  ±0.025  mm  and  ±0.076  mm,  respectively.  The  velocity 
measurements  were  considered  accurate  to  within  ±3  m/s.  The 
accuracy  imposed  on  the  velocity  measurement  system  is  taken 
from  an  analysis  of  the  uncertainty  of  data  generated  by  the 
pressure  transducers.  It  was  assumed  for  this  analysis  that  the 
wave  forms  recorded  provided  the  optimum  level  of  resolution  and 
the  information  was  collaborated  by  the  infrared  detector  and 
the  high  speed  photography  techniques  previously  discussed. 
Correlating  the  velocity  data  between  the  various  measurement 
systems  provided  an  increased  level  of  confidence  in  the 
accuracy  of  the  measurement.  By  choosing  the  results  of  a 
particular  impact  test  as  a  baseline,  the  effect  of  possible 
inaccuracies  in  the  experimental  parameters  on  the  computed 
yield  stress  can  be  studied. 

Table  5  shows  a  set  of  baseline  values  for  two  different 
pure  copper  specimens  with  aspect  ratios  of  1.5  and  10.  By 
examining  the  data  from  specimens  with  different  geometries,  the 
relative  influence  of  possible  inaccuracies  can  be  considered. 
Perturbed  values  of  these  parameters,  according  to  the  assumed 
accuracy  of  the  measurement,  are  shown  in  the  columns  headed  with 
a  plus  and  a  minus  sign. 

In  Table  6,  the  dynamic  yield  stress  values  are  given  for 
the  baseline  and  from  computations  using  the  perturbed  parameter 
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TABLE  5. 


EXPERIMENTAL  TEST  DATA 


ASPECT  RATIO 

T  .5 

TO 

MEASUREMENT 

- 

BASELINE 

+ 

- 

BASELINE 

Lf 

9.093  mm 

9.119  mm 

9.144  mm 

60.528  mm 

If 

3.861  mm 

27.000  mm 

27.076  mm 

27.153  mm 

V 

165.3  m/s 

168.3  m/s 

171.3  m/s 

163.0  m/s 

166.0  m/s 

169.0  m/s 

values.  In  making  the  calculations  only  one  of  the  three 
measurements  were  varied  from  the  baseline  value.  For  example, 
in  the  first  row,  using  a  specimen  with  L/D  of  1.5,  the  yield 
stress  was  calculated  using  a  value  of  the  final  length  of 
9.093  mm,  9.119  mm,  and  9.144  mm;  the  baseline  values  were  used 
for  the  other  two  parameters,  if  and  V.  In  the  second  and  third 
row  of  Table  6,  the  perturbed  parameter  was  the  undeformed  length 
and  velocity,  respectively.  A  quantitative  estimate  of  the 
uncertainty  can  be  determined  from  a  ratio  of  the  difference 
between  the  computed  yield  stresses  found  using  the  perturbed 
experimental  parameters  to  the  baseline  value.  The  percent 
errors  were  3.4,  2.0,  and  7.2,  for  measurements  of  Lf,  if,  V, 
respectively,  for  a  specimen  with  an  aspect  ratio  of  1.5.  By 
comparison,  the  second  set  of  data  produced  percentage  error 
values  of  0.6,  0.4,  and  7.2,  from  measurements  of  Lf,  if,  and  v, 
respectively.  The  velocity  measurement  provided  the  greatest 
source  of  uncertainty  in  both  data  sets.  By  comparison  to  the 
velocity  measurement,  the  deformed  specimen  geometry  had  little 
influence  on  the  overall  uncertainty  of  the  computed  yield  stress 
value . 

An  argument  could  be  made  that  additional  combinations  of 
the  various  parameters  would  produce  a  larger  deviation  from  the 
baseline  value.  Note  in  Table  6  that  a  reduction  of  the  geometry 
parameters,  Lf  and  if,  had  an  opposite  effect  on  the  computed 
yield  stress  value.  Reducing  the  final  undeformed  length,  If, 


increases  the  stress  value,  while  a  reduction  in  the  final 
length,  Lf,  decreases  the  calculated  stress.  It  was  considered 
that  the  actual  uncertainties  in  postmortem  measurements  would 
be  small,  or  possibly  cancel  one  another,  by  comparison  to  those 
found  in  the  velocity  measurement.  Therefore,  the  yield  stress 
values  reported  in  Table  6  as  baseline  values  and  those  given  in 
Appendix  B,  are  accurate  to  within  ±5  percent. 


TABLE  6.  COMPUTED  YIELD  STRENGTH  VALUES 


ASPECT  RATIO 

1  .5 

io 

VARIANCE 

- 

BASELINE 

+ 

- 

BASELINE 

4- 

L, 

549  Mpa 

559  Mpa 

568  Mpa 

510  Mpa 

512  Mpa 

513  Mpa 

If 

564  Mpa 

559  Mpa 

553  Mpa 

513  Mpa 

512  Mpa 

51 1  Mpa 

V 

539  Mpa 

559  Mpa 

579  Mpa 

493  Mpa 

512  Mpa 

530  Mpa 

5.  EFFECT  OF  SPECIMEN  GEOMETRY 

The  influence  of  initial  specimen  geometry  on  the  calculated 
yield  stress  was  another  effect  of  interest  in  this 
investigation.  Figure  15  shows  the  calculated  value  of  yield 
stress  as  a  function  of  the  aspect  ratio,  L/D,  of  the  specimen. 
These  data  show  that  the  specimen  geometry  did  not  influence  the 
yield  stress  value  calculated  by  the  model,  with  the  exception 
being  the  2024-T4  material. 

The  2024-T4  aluminum  shows  an  increase  in  the  yield  stress 
with  short  specimens.  The  value  of  the  yield  stress  varies  from 
750  MPa,  for  L/D’s  of  5  and  greater,  to  near  900  MPa,  for  L/D's 
of  1.5  and  2.  While  the  uncertainty  is  somewhat  greater  for 
smaller  specimens,  in  terms  of  the  yield  calculations,  there  are 
sufficient  data  to  suggest  the  indication  of  a  specimen  geometry 
related  trend. 


43 


< 

ct : 


cj 

UJ 

r\ 

in 

< 

tn 

f—  > 
21 

UJ  H 

ii 

£l5 

CL  (V 
X  p 
UJOO 

& 


o 

s< 

<  U 

Q 


2  2 
2  z) 
—  — 

Of  ■'■'  5  ID 

~CL<  < 
CL  Q_ 

O  (j  H*  I— 

I  I 

O  UJ  'T  r~ 

I  I —  CN  CO 
Lu  Cl  o  o 
O  Q  CN  CD 


0*0*03  + 


O 

o 

o 

o 

o 


♦  40+ 


o  o 


00  o 

o  o 


a 

o#  *a  + 


♦  oa**+ 

CJ  0  ♦ 


-i — 1 — | — i — r 


1  1  I 


PS 

P 

P 


L  O 

ho 

CO 


o  <r 

0 ,  v 

C£>  1 - 


o 

o  UJ 
CL 
00 
< 


b  N- 


o 

■o 

b  CN 


0 

O 

O 

O 

O 

O 

O  O 

0 

O 

O 

O 

O 

O 

CN 

O 

OO 

IO 

T 

CN 

Ddn)  ‘H10N3yiS  G13IA  OllAVNAQ 


44 


Figure  15.  Values  of  the  Yield  Stress  as  a  Function  of  the 
Specimen  Aspect  Ratio 


If  the  increase  in  yield  stress  occurred  because  of  inaccura¬ 
cies  in  the  postmortem  measurements,  the  results  for  the  other 
three  materials  should  be  affected  as  well.  This  supports  the 
indication  that  the  specimen  geometry  does  influence  the  yield 
stress  for  this  type  of  material.  How  the  specimen  geometry 
affected  the  results  can  be  studied  by  examining  the  data 
obtained  from  the  experiment. 

Recall  that  the  nondimensional  forms  of  the  final  length  and 
final  undeformed  length  were  shown  to  decrease  linearly  with 
increasing  velocity,  Figure  14.  The  statistical  analysis,  Table 
3,  shows  that  the  curve  fit  by  linear  regression  of  the  data 
produced  a  slope  that  was  equal  for  the  two  parameters,  with  the 
exception  being  the  2024-T4  alloy.  Data  for  the  2024-T4 
material  show  that  the  final  undeformed  length  changed  at  a 
faster  rate,  with  velocity,  than  did  the  final  length  of  the 
specimen.  These  results  indicate  direct  physical  evidence  of  a 
unique  material  behavior  for  this  alloy. 

One  of  the  methods  used  to  characterize  material  behavior  by 
Whiff in,  Reference  9,  in  the  original  investigation  was  to  plot 
the  nondimensional  form  of  the  final  undeformed  length  against 
the  nondimensional  form  of  the  final  length,  Figure  16, 

Reference  9.  For  the  current  investigation  similar  data  are 
plotted  in  Figures  17(a)  and  17(b).  Whiffin  used  this  method  of 
displaying  the  data  to  comply  with  a  requirement  of  the  Taylor 
model,  from  which  the  value  derived  from  an  explicit  formula  of 
the  yield  stress,  Equation  (18) ,  page  10,  could  be  used  to  calcu¬ 
late  a  more  precise  value  using  the  correction  contours. 

The  purpose  for  showing  Figures  16  and  17,  is  to  discuss  the 
results  of  the  original  investigation  in  light  of  the  current 
work.  The  results  shown  in  Figure  16  for  Normal  Dural  are  to  be 
compared  with  the  results  shown  in  Figure  17  for  the  2024-T4 
material . 
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The  velocities  shown  next  to  the  data  points  in  Figure  16 
indicate  that  the  complex  state  of  stress  in  the  material  accom¬ 
panying  fracture  was  not  considered  in  that  investigation. 
Therefore,  only  one  of  the  four  data  points  shown  for  Normal 
Dural  would  be  acceptable  in  regard  to  the  no  fracture  require¬ 
ment  imposed  on  the  current  test  data.  The  data  point  for  impact 
velocity  of  640  ft/s,  (195  m/s),  is  in  agreement  with  the  values 
in  the  current  investigation. 

Whiffin  draws  the  conclusion  from  the  data  that  for  a 
specific  value  of  lf/L,  the  ductility  of  various  materials  can 
be  identified  from  the  parameter  Lf/L.  Whiffin  identifies  the 
trend  that  at  a  specific  value  of  lf/L  higher  values  of  Lf/L 
occur  for  more  ductile  materials.  In  Figure  17,  it  can  be  seen 
that  for  a  specific  value  of  the  final  undeformed  length 
parameter,  lf/L,  increased  ductility  will  decrease  the  value  of 
the  final  length,  not  increase  it.  For  example,  the  2024-T4 
material  has  less  ductility  than  does  6061-T6,  and  is, 
therefore,  found  to  the  right  of  the  data  for  6061-T6.  This 
misinterpretation  of  the  data  occurred  because  of  the  limited 
amount  of  data  that  was  available  and  the  relatively  crude 
measurement  technique  employed  to  make  the  postmortem 
measurements . 

A  trend  in  the  data  can  be  identified,  for  the  2024-T4 
alloy,  when  the  aspect  ratio  of  the  specimen  is  plotted  on  the 
data  point,  Figure  17(b).  In  general,  the  lower  L/D's  are 
shifted  toward  higher  values  of  the  nondimensional  final  length 
parameter.  This  trend  in  the  data  identifies  the  relationship 
between  the  postmortem  geometry  measurements  that  are  influencing 
the  computed  value  of  the  yield  stress.  As  shown  by  the  sensiti¬ 
vity  analysis,  reducing  the  final  undeformed  length  produces  an 
increase  in  the  yield  stress.  Consequently,  this  shift  in  the 
data  for  specimens  with  lower  L/D's  produced  the  results  shown  in 
Figure  15. 
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MATERIAL 


REF. 

NO. 


STATIC  STRENGTH  AVERAGE  DYNAMIC 
(TONS/SQ.  IN.)  STRENGTH  (TONS/SQ.  IN.) 


1 

ARMOUR  PLATE 

73 

115 

2 

•NORMAL*  DURAL 

19 

31 

3 

HARD  DURAL 

25-85 

39 

4 

MILE  STEEL 

18 

49 

5 

ARM  CO  IRON 

— 

50 

6 

SOR  DURAL 

8-15 

19 

7 

LEAD 

— 

2.8 

8 

COPPER 

— 

15.5 

.  THE  SMALL  FIGURES  IN  THE  DIAGRAM  SHOW  THE  STRIKING  VELOCITY  EMPLOYED. 


Figure  16.  Nondimensional  Form  of  the  Undeformed  Final 
Length  as  a  Function  of  the  Nondimensional 
Final  Length  Parameter 


U/  L 


^^-^ure  17a.  Ballistic  Geometry  Data,  in  Nondimensional 
Form,  for  the  Current  Investigation 
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TAYLOR  IMPACT  EXPERIMENT 
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Figure  17b.  Ballistic  Geometry  Data,  in  Nondimens ional 
Form,  for  2024-T4  Specimens 
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6 .  MUSHROOM  GROWTH 

The  measurement  of  the  final  diameter  of  the  specimen  at  the 
target/ specimen  interface  indicates  that  the  mushroom  growth  is 
correlated  with  the  available  energy  on  impact,  Figure  18.  A 
comparison  between  short  and  long  specimens  shows  that  the  final 
diameter  is  proportionally  greater  for  longer  specimens, 
depending  on  the  mass  increase,  for  equal  impact  velocities. 
These  results  indicate  that  the  final  diameter  of  the  mushroom 
is  dependent  on  the  strength  of  the  material  and  the  kinetic 
energy  at  impact. 
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I AYLOR  IMPACT  EXPERIMENT 
0.60  FNAL  DIAMETER  vs  IMPACT  VEIOCIIT 
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Figure  18.  Final  Diameter  as  a  Function  of  the  Impact  Velocity 


7. 


HIGH  SPEED  PHOTOGRAPHY 


The  Cordon  framing  camera  was  used  to  generate  a  photographic 
record  of  a  particular  impact  test  (UK-145) .  The  camera  was 
operated  at  a  framing  rate  of  300,000  frames  per  second.  At  this 
rate,  the  resolution  of  the  photographic  data  was  3.3  3  /is.  The 
impact  event  was  photographed  using  a  shadowgraph  technique. 
Measurements  taken  directly  from  the  photographs  provided  data  on 
the  mushroom  growth  and  the  rear  end  position  both  as  a  function 
of  time. 

The  data  obtained  on  the  mushroom  growth  rate,  Figure  19, 
indicates  the  magnitude  of  the  strain  rates  seen  at  the 
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Figure  19.  Strain  at  the  Target/Specimen  Interface 
as  a  Function  of  Time 
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specimen/target  interface.  The  photographs  show  that  the  bulk  of 
the  mushroom  growth  process  occurs  within  the  first  few  microse¬ 
conds  and  is  completed  by-  40  ns  after  impact.  The  slope  of  the 
curve  in  Figure  19  gives  a  strain  rate  of  7.5  x  104/s  during  the 
first  7  ms. 

The  data  on  the  position  of  the  rear  end  of  the  specimen  as 
a  function  of  time  shows  that  it  took  120  ms  to  bring  the 
specimen  to  rest,  Figure  20.  The  velocity  of  the  rear  end  is 
simply  the  slope  of  the  curve.  The  data  show  that  the  initial 
velocity,  190  m/s,  did  not  change  until  approximately  45  jis 
after  impact. 


Figure  20.  Position  of  the  Undeformed  End  of  the 
Specimen  as  a  Function  of  Time 
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The  analytical  models  which  attempt  to  describe  the  physical 
process  assume  that  the  plastic  front  moves  at  a  uniform  velocity 
throughout  the  event.  The  calculated  plastic  wave  speed  from  the 
a/0  model  should  be  considered  to  be  an  average  value.  For  the 
particular  test  that  was  photographed  with  the  framing  camera,  a 
comparison  of  the  plastic  wave  speed  calculated  from  the  model 
can  be  made  with  experimental  results.  Knowing  the  duration  of 
the  event  and  the  distance  traveled  by  the  plastic  wave,  a  velo¬ 
city  was  determined.  For  test  UK-145,  the  model  computed  a  value 
of  202  m/s  for  the  parameter,  \.  From  the  photographs,  the  plas¬ 
tic  wave  speed  is  calculated  to  be  212  m/s,  which  is  in  good 
agreement  with  the  value  predicted  by  the  model. 
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SECTION  V 


MICROSTRUCTURE 


1 .  INTRODUCTION 

The  materials  considered  for  this  investigation  varied 
widely  in  their  microstructural  features.  The  differences  in 
deformed  geometries  reflect  how  the  microstructure  of  each 
material  influenced  the  mechanisms  for  plastic  deformation.  The 
principal  aspect  of  this  study  is  how  the  basic  microstructure 
affects  the  way  the  material  absorbs  the  energy  of  impact. 

Macroscopically ,  the  features  of  a  deformed  specimen  were 
similar  for  all  of  the  materials.  The  grains  on  the  axis  near 
the  anvil,  which  are  subjected  to  high  compressive  loads,  are 
extensively  deformed.  The  longitudinal  compressive  strains 
impart  radial  motion  to  the  deforming  body  to  produce  the 
characteristic  mushroomed  profile.  The  extent  of  the  deforma¬ 
tion  process  is  governed  by  various  aspects  of  the  microstruc¬ 
ture,  which  is  controlled  by  the  thermo-mechanical  history  of 
the  material  and  its  chemical  composition. 

The  structures  found  in  the  materials  used  for  this 
investigation  were  of  three  types:  single  phase,  multiphase 
containing  inclusion  particles,  and  multiphase  containing 
inclusions  and  coherent  precipitates.  The  first  of  these 
structures  is  found  in  OFHC  copper  (99.99  Cu) .  The  second  type 
is  found  in  DPTE  copper  and  in  6061-T6  aluminum.  The  last  type 
is  found  in  2024-T4  aluminum.  The  constituents  found  in  each  of 
the  materials  are  given  in  Table  7.  The  heat  treatments 
specified  for  the  aluminum  materials  correspond  to  a  solution 
anneal  followed  by  artificial  aging,  for  6061,  and  natural 
aging,  for  2024.  Some  specimens  of  OFHC  copper  were  annealed 
for  1  hour  at  600°C  in  a  vacuum.  The  nomenclature  annealed 
copper  and  half  hard  copper  will  henceforth  be  used  to  denote 
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the  CFHC  copper  in  the  annealed  and  as-received  conditions, 
respectively. 

All  of  the  materials  were  examined  by  optical  microscopy 
both  before  and  after  impact.  The  OFHC  coppers  were  more 
extensively  studied.  In  addition  to  optical  examination,  X-ray 
diffraction  and  transmission  electron  microscopy  (TEM) 
techniques  were  used  to  study  the  microstructure  and  the  effect 
of  high  strain  rate  deformation. 


DPTE  COPPER 


2  . 

The  DPTE  copper  material  is  a  copper/tellurium  alloy.  Tellu¬ 
rium  is  added  as  a  solution  strengthening  element  in  small 
amounts.  A  small  concentration  of  phosphorous  is  added  as  a 
deoxidizer.  The  phosphorous  removes  oxygen  from  the  copper 
matrix  by  forming  oxide  inclusion  particles.  During  hot  rolling 
of  the  material  these  inclusions  are  broken  up  and  elongated  in 
the  rolling  direction.  The  final  structure  after  cold  working 
the  material  shows  a  fine  grain  structure  with  an  extensive 
number  of  stringers  elongated  parallel  to  the  rod  axis,  Figure 
21. 


Figure  21.  Optical  Micrograph  Showing  the  Microstructure 
of  DPTE  Copper 

3.  6061-T6  ALUMINUM 

The  6061-T6  is  a  ternary  aluminum  alloy  containing  small 
concentrations  of  silicon  and  magnesium.  Age  hardening  will 
produce  a  matrix  of  aluminum  with  magnesium  in  solution  and 
precipitates  of  Mg2Si.  These  precipitate  particles  will  be 
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spherical  in  shape.  Figure  22,  and  will  be  incoherent  with  the 
matrix.  In  commercial  grades,  inclusions  of  aluminum  oxides  and 
iron/silicon  compounds  will  be  found  that  are  not  affected  by 
subsequent  heat  treatment.  The  size  of  the  misfit  between 
lattice  types,  the  aluminum  matrix  and  the  precipitate,  produces 
a  small  internal  strain  field  surrounding  each  precipitate.  This 
small  strain  field  produces  little  disruption  to  the  motion  of 
dislocation  which  occurs  under  an  applied  load. 


Figure  22.  Optical  Mici'ograph  Showing  the  Microstructure 
of  6061-T6  Aluminum 

The  shape  of  the  precipitate  particles  reduces  the  number  of 
potential  sites  for  localized  stress  concentration.  In  turn, 
preventing  the  localized  buildup  of  stress  inhibits  crack 
nucleation  and  growth.  This  microstructure  in  the  6061-T6  alloy 
provides  a  relatively  low  yield  strength,  but  gives  it  a  high 
toughness,  i.e.,  the  ability  to  absorb  energy  before  fracture. 
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4.  2024-T4  ALUMINUM 

The  2024-T4  is  a  ternary  aluminum  alloy  containing  copper 
and  magnesium.  Allowing  the  material  to  naturally  age  after 
solution  anneal  provides  the  opportunity  for  coherent  Al2Cu 
precipitates  to  form  on  (001)  habit  planes.  In  addition, 
magnesium  and  aluminum  oxides  form  inclusion  particles  in  the 
matrix,  Figure  23.  The  formation  of  the  coherent  Al2Cu 
precipitates,  which  have  a  large  misfit  with  the  matrix, 
prevents  the  easy  motion  of  dislocations  on  the  slip  planes. 
Therefore,  the  material  shows  increased  strength  by  comparison 
to  the  6061-T6  alloy,  see  Table  1. 


Figure  23.  Optical  Micrograph  Showing  the  Microstructure 
of  2024-T4  Aluminum 

5.  OFHC  COPPER 

Optical  examination  of  the  half  hard  OFHC  copper  material 
shows  that  slip  is  the  primary  mode  of  accommodating  the 
deformation  introduced  by  drawing,  Figure  24.  TEM  thin  foils  of 
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the  material  showed  an  extensive  substructure  containing  an 
inhomogeneous  distribution  of  dislocations,  Figure  25.  The 
parallel  bands  seen  in  Figure  25  indicate  that  deformation  twins 
were  formed  during  the  processing  of  the  material. 


Figure  24.  Optical  Micrograph  of  the  Microstructure  of 
Half  Hard  OFHC  Copper 

Reports  by  Ahlborn  and  Wassermann,  Reference  11,  and 
Wassermann,  Reference  12,  cite  evidence  that  the  following 
conditions  give  an  increased  incidence  of  mechanical  twinning  in 
drawn  copper  wire: 

a.  lowering  the  deformation  temperature,  Reference  13, 

b.  increasing  the  amount  of  deformation, 

c.  lowering  the  stacking  fault  energy  by  alloying. 

The  as  received  material  was  in  a  half-hard  condition,  which 
means  that  following  the  final  anneal  there  was  a  reduction  in 
cross  sectional  area  of  30  to  60  percent.  This  amount  of 
straining  satisfies  condition  b,  reported  above,  and  is 
confirmed  by  the  TEM  observations  in  Figure  25. 
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Pole  figures  were  generated  to  determine  the  textures  of  the 
half-hard  and  annealed  OFHC  copper  materials.  Because  of  the 
drawing  operation,  the  texture  in  the  half  hard  material  was 
found  to  be  oriented  parallel  to  the  rod  axis  and  was  a 
superposition  of  (111)  and  (001)  fiber  textures.  The  annealed 
material  had  a  random  grain  orientation  with  an  approximate 
grain  size  of  40  ^m.  Compare  Figures  24  (half  hard)  and  26 
(annealed) . 

6.  MICROSTRUCTURE  RESULTING  FROM  IMPACT 

With  different  metallurgical  histories,  two  copper  specimens 
which  impact  the  target  at  similar  velocities  accommodated  the 
energy  in  significantly  different  ways.  Examination  of  the 
post-impact  geometries  reveal  how  the  microstructure  influenced 
the  deformation,  Figure  27.  The  specimen  on  the  left  was 
deformed  in  the  half-hard  condition.  The  specimen  on  the 
right  was  deformed  in  the  annealed  condition.  The  three  post¬ 
mortem  measurements  of  the  geometry  that  were  visibly  influenced 
by  the  microstructural  differences  are  the  final  mushroom 
diameter,  the  final  length,  and  the  extent  of  the  plastic 
deformation  in  the  specimen. 

The  half-hard  copper  specimen  had  a  much  greater  change  in 
diameter.  The  annealed  specimen  had  a  greater  reduction  in  final 
length  and  was  traversed  entirely  by  the  plastic  wave  in  the 
material.  From  these  observations,  it  is  quite  evident  that  the 
dislocation  density  in  the  material  had  a  significant  role  in 
how  the  impact  energy  was  absorbed. 

Following  impact,  both  types  of  copper  specimens  were 
sectioned  parallel  to  the  rod  axis.  An  optical  examination  was 
conducted  near  the  mushroom  face  close  to  the  original  rod  axis, 
Figures  28  and  29.  The  large  compressive  loads  are  evident  by 
the  collapse  of  the  grains.  After  impact,  the  grains  of  the 
annealed  material  run  parallel  to  the  target  face,  Figure  28. 
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Figure  27.  Impact  Specimens  of  OFHC  Copper 

In  contrast,  the  half-hard  copper  specimen  contains  small, 
randomly  oriented  grains,  Figure  29.  These  grains  indicate  that 
the  mechanical  history  of  the  material  created  a  condition 
suitable  for  recrystallization  to  occur  as  a  result  of  impact. 

To  confirm  this  observation,  back  reflected  Laue  patterns 
were  generated  from  both  types  of  impact  specimens.  The  x-ray 
beam  was  positioned  to  strike  near  the  same  location  as  shown  in 
the  optical  micrographs.  Results  for  the  annealed  impact 
specimens  showed  two  diffuse  rings.  The  impact  specimen  used  in 
the  half-hard  condition  showed  two  rings  made  up  of  numerous 
intense  spots.  These  bright  spots  occur  because  of  the  increase 


Figure  29. 


0.4  mm 

Optical  Micrograph  Showing  the  Microstructure 
of  Half  Hard  OFHC  Copper  Specimen  After  Impact 
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in  size  of  the  coherently  diffracting  domain  created  by  the 
recrystallization  process. 

The  high  dislocation  density  in  the  half  hard  material  had  a 
significant  effect  on  the  behavior  in  both  impact  tests  and 
tensile  tests.  In  tensile  tests,  performed  for  quasi-static 
characterization  of  the  material,  the  high  dislocation  density 
effectively  eliminated  the  strain  hardening  region  of  the  load 
versus  time  plot,  Appendix  D.  The  load  increased  linearly  until 
localized  plastic  deformation  occurred  so  that  the  yield  and 
ultimate  strengths  of  the  material  coalesced  into  one  point. 
During  the  deformation  on  impact,  the  dislocation  density 
influenced  how  the  material  absorbed  the  energy.  As  previously 
mentioned,  this  observation  is  evident  by  the  macroscopic 
features  of  the  post-impact  geometry,  Figure  27.  The 
post-impact  geometries  suggests  that  recrystallization  occurred 
in  the  half-hard  material  because  the  heat  generated  by  plastic 
work  was  more  intense.  In  addition,  the  increase  in  stored 
energy  in  the  material  caused  by  cold  working  reduces  the  amount 
of  energy  required  for  recrystallization.  The  temperature  for 
the  recrystallization  of  copper  is  225°C. 

After  impact,  TEM  thin  foils  were  made  from  both  types  of 
copper  materials.  The  half-hard  material  showed  a  high 
incidence  in  the  number  of  deformation  twins,  Figure  30.  Since 
the  material  had  contained  twins  initially,  the  source  of  the 
twins  could  have  been  either  the  impact  event  or  the  processing 
of  the  original  rod  stock.  However,  deformation  twins  were  also 
found  in  the  specimens  of  annealed  material.  These  could  only 
have  been  produced  by  impact. 

TEM  replicas  were  taken  from  the  surface  of  the  sectioned 
half-hard  impact  specimen.  Examination  of  these  replicas 
revealed  an  orientation  relationship  that  identified  the 
formation  of  a  twin  with  the  impact  event,  Figure  31.  Previous 
investigators,  Reference  14,  identified  similar  artifacts  as 
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compression  twins  formed  in  shock  loaded  single  crystals  of  pure 
copper.  (See  Figure  14,  in  Reference  14.) 


Figure  30.  TEM  Micrograph  Showing  Deformation  Twins  in  OFHC 
Copper 

• 

Brilhart  and  coworkers.  Reference  15,  investigated  the 
formation  of  deformation  twins  in  polycrystalline  pure  copper  by 
shock  loading  the  material  using  flyer  plate  experiments.  In 
these  experiments  the  pressures  responsible  for  the  formation  of 
deformation  twins  were  on  the  order  of  7.5  GPa. 

The  observation  of  deformation  twins  in  pure  copper  Taylor 
test  specimens  suggests  the  level  of  pressure  in  the  material  on 
impact.  The  pressure  on  contact  between  the  specimen  and  the 
target  was  est: mated  by  a  model  constructed  using  Hugoniot 
relationships,  as  detailed  in  Appendix  E.  From  this  analysis  the 
pressure  at  the  moment  of  impact,  in  a  copper  specimen  striking 
a  steel  anvil  at  an  initial  velocity  of  200  m/s,  is  3.9  GPa. 

This  estimated  pressure  is  nearly  an  order  of  magnitude  above 
the  dynamic  yield  stress  of  the  material. 
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Figure  31.  Photo-Micrograph  of  a  TEM  Replica  Showing  a 

Shock  Formed  Twin  in  Half-Hard  OFHC  Copper  After 
Impact 

The  mechanical  behavior  of  copper  was  significantly  affected 
by  increasing  the  impurity  content.  Increasing  impurity  content 
reduced  ductility.  In  static  tests,  the  percent  elongation,  or 
strain,  at  which  failure  occurred  was  lower  in  the  DPTE  copper 
than  in  the  OFHC  copper. 

In  dynamic  tests,  the  data  showed  that  the  DPTE  copper 
failed,  by  fracture  on  the  periphery  of  the  mushroom  diameter,  at 
impact  velocities  greater  than  160  m/s.  Failure  of  the  impact 
specimen  resulted  from  the  large  circumferential  strains  produced 
by  the  radial  motion  of  the  deforming  material.  Among  DPTE 
copper  specimens  which  had  not  cracked  radially,  the  largest 
mushroom  diameter  measured  was  20  percent  less  than  for  OFHC 
copper  specimens. 

The  type  of  fracture  found  in  the  DPTE  copper  material  was 
indicative  of  the  microstructure.  The  high  level  of  impurities 
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provided  numerous  locations  within  the  material  for  the  localized 
build  up  of  stress.  The  ability  of  t'  e  matrix  to  accommodate  the 
strain  was  reduced  by  the  addition  of  tellurium.  Numerous  sites 
for  stress  concentrations,  and  reduced  matrix  ductility,  produced 
a  material  that  fractured  at  relatively  low  strain  levels.  The 
mode  of  failure  was  by  brittle  fracture  as  indicated  by  the 
grainy  texture  of  the  fracture  surfaces. 

By  comparison,  the  OFHC  copper  specimens  failed  in  a  ductile 
manner.  Close  observation  of  the  mushroom  periphery  in  Figure 
27  shows  the  early  stages  of  crack  growth.  Note  the  region  of 
deformation  which  surrounds  each  crack  and  is  typical  of  ductile 
failure . 

The  inclusions,  or  stringers,  in  the  DPTE  copper  material 
provide  a  unique  way  of  examining  the  results  of  deformation 
within  the  specimen.  During  the  original  processing  of  the  rod 
stock  the  oxide  inclusions  were  broken  and  elongated  in  a 
direction  parallel  with  the  axis  of  the  rod.  After  impact, 
examination  of  the  new  orientation  of  the  stringers  and  grain 
boundaries  within  the  material  revealed  certain  aspects  about  the 
deformation  resulting  from  impact.  During  impact,  the  grains 
compress  in  the  direction  of  the  specimen  axis,  while  expanding 
radially.  Comparison  of  the  grain  boundaries  before  and  after 
impact.  Figures  21  and  32,  reveals  that  the  grains  at  the  axis 
now  run  parallel  to  the  anvil.  To  accommodate  the  large 
compressive  strains,  the  free  boundary  moves  radially  outward. 

In  the  region  nearest  the  free  surface  stringers  that  are 
initially  oriented  axially  reorient  as  the  free  surface  changes 
shape,  Figure  33.  Therefore,  the  shortening  of  these  stringers 
is  less  than  those  near  the  specimen  axis,  allowing  these 
inclusion  stringers  to  retain  their  basic  shape.  The  radial 
expansion  that  has  occurred  in  the  material  is  depicted  by  the 
new  alignment  of  the  stringers.  Compare  Figures  21,  32,  and  33. 
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The  two  aluminum  materials  responded  to  the  deformation  in  a 
manner  consistent  with  their  respective  microstructures.  The 
grains  on  the  axis  of  the  6061-T5  aluminum  alloy  readily- 
deformed,  Figure  34.  However,  the  2024-T4  material,  Figure  35, 
resisted  the  large  amount  of  deformation  seen  in  the  6061-T6 
alloy.  This  is  evidence  of  what  effect  second  phase  particles 
have  on  the  dislocations  as  opposed  to  incoherent  intermetallic 
particles.  The  finely  dispersed  second  phase  precipitates  in 
the  2024-T4  alloy  act  as  pegs  for  the  dislocations,  which 
attempt  to  move  along  the  slip  planes  in  response  to  the  applied 
stress.  It  can  be  seen  in  Figure  35  that  slip  bands  have  been 
formed  in  certain  grains.  These  grains  have  crystallographic 
orientation  which  maximizes  the  shear  stress  found  on  the  slip 
plane.  Once  the  Al2Cu  precipate  is  sheared,  additional 


Figure  34.  Optical  Micrograph  Showing  the  Microstructure 
in  6061-T6  Aluminum  After  Impact 

dislocations  readily  move  through  the  material.  The  inclusion 
particles  of  6061-T6  aluminum  being  larger  in  size  and  spaced 
further  apart  than  the  precipitates  in  2024-T4  provides  little 
resistance  to  the  motion  of  dislocations. 
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The  energy  which  can  be  absorbed  by  a  material  is  defined  by 
its  toughness  and  is  indicated  in  the  Taylor  anvil  test  by  the 
strain  to  failure  and  the  extent  of  the  plastic  deformation. 

The  2024-T4  alloy  has  high  strength,  but  relatively  low  tough¬ 
ness,  compared  with  the  6061-T6  alloy.  The  geometry  of  the 
inclusion  particles  found  in  2024-T4  are  very  irregular,  Figure 
23.  This  shape  produces  localized  stress  concentrations  which 


Figure  35.  Optical  Micrograph  Showing  the  Microstructure 
in  2024-T4  Aluminur  After  Impact 

lead  to  crack  nucleation.  Once  initiated,  the  crack  propagates 
through  the  material  at  the  interface  between  the  inclusion 
particles  and  the  matrix.  In  contrast,  the  6061-T6  alloy  is 
relatively  tough  because  of  the  spherically  shaped  particles, 
which  c’o  not  s ignif icantly  increase  the  local  stress  level, 
Figure  22.  Similarly,  the  DPTE  copper  had  low  toughness,  while 
the  OFHC  copper  could  absorb  an  extensive  amount  of  energy  by 
plastic  deformction. 
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SECTION  VI 


TWO-DIMENSIONAL  MODELING 


1.  INTRODUCTION 

One-dimensional  interpretations  of  the  Taylor  test  have  been 
discussed  previously.  In  addition,  the  Taylor  test  can  provide 
experimental  information  to  be  used  with  other  types  of  computer 
models.  These  models  consist  of  two-dimensional  simulations  of 
the  impact  event.  The  key  element  in  such  computer  programs  is 
the  constitutive  relationship  that  predicts  the  material  response 
in  the  event. 

2.  CONSTITUTIVE  RELATIONSHIPS 

The  constitutive  relationship  needed  in  the  model  must 
provide  a  relationship  between  the  independent  variables,  i.e., 
strain,  strain-rate,  temperature,  and  the  dependent  variable, 
stress.  Often,  these  terms  are  combined  in  a  generic  formula 
along  with  various  material  dependent  adjustment  factors.  The 
adjustment  factors  will  take  into  consideration  the  sensitivity 
of  the  stress  level  in  the  material  under  various  conditions. 

For  example,  a  common  term  used  in  this  type  of  characterization 
is  work  hardening  of  the  material.  Typically,  this  strain 
dependent  term  has  a  coefficient,  which  dictates  the  amount  of 
work  hardening,  and  an  exponent  term,  which  governs  the  rate  of 
work  hardening.  The  coefficient  and  the  exponent  are  material 
dependent  factors  within  the  equation.  One  example  of  such  an 
expression,  used  with  some  success  by  investigators  Johnson  and 
Cook,  Reference  4,  is  as  follows: 

a  =  (A  +  Ben)  (1  +  Cline)  (1  -  T*m)  (42) 

where  A,  B,  C,  n,  m,  are  material  parameters  which  must  be 
determined  in  laboratory  test  for  the  material  under 
investigation.  The  parameter  T*  is  the  homologous  temperature. 
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The  first  term  on  the  right-hand  side  of  Equation  (42)  is 
the  expression  given  by  Ludwik,  Reference  16,  for  describing  the 
parabolic  work  hardening  behavior  of  a  metal.  The  second  and 
third  terms  contain  the  influence  of  strain  rate  and  tempera¬ 
ture,  respectively,  on  the  stress-strain  curve. 

Other  investigators,  Zerilli  and  Armstrong,  Reference  5, 
have  modified  the  equation  proposed  by  Johnson  and  Cook.  The 
investigation  by  Zerilli  and  Armstrong  provided  an  enhanced 
material  behavior  description  for  use  in  multi-dimensional 
models.  It  was  pointed  out  that  Equation  (42)  was  an 
improvement  over  previously  used  constitutive  relationships, 
e.g.,  elastic/plastic,  or,  elastic/plastic  with  linear  strain 
hardening.  However,  Equation  (42)  was  deficient  in  describing 
certain  aspects  of  observed  material  behavior.  Factors  known  to 
have  an  affect  on  the  material  behavior,  i.e.,  grain  size, 
dislocation  density,  etc.,  were  noticeably  absent.  An  equation 
proposed  by  Zerilli  and  Armstrong  for  face  centered  cubic  (fee) 
material  is  as  follows: 

a  =  &OQ  +  kJi-0  •  5  +  C2  €  0  •  5exp  (-Cj  T  +  C4  T!l.ne  )  (43) 

where  k  is  the  microstress  intensity  at  the  grain  boundary  and  1 
is  the  average  grain  size.  The  term  a aG  is  a  stress  component 
associated  with  the  effect  of  solute  atoms  and  the  original 
dislocation  density  in  the  material.  The  influence  of  the 
microstructure  is  incorporated  into  the  first  and  second  terms 
on  the  right  hand  side  of  Equation  (43) ,  as  well  as  in  the 
material  dependent  coefficients,  C2  ,  C3  ,  C4  .  The  last  term 
combines  the  influences  of  strain,  strain  rate,  and  temperature 
on  the  material  behavior. 

A  similar  expression  was  constructed  for  body  centered  cubic 
(bcc)  materials.  This  expression  was  modified  slightly  from 
Equation  (43)  by  separating  the  strain  hardening  term  from  the 
strain  rate  and  temperature  components.  As  suggested  fc>y  Zerilli 
and  Armstrong,  Reference  5,  this  modification  is  in  agreement 
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with  observed  behavior  of  bcc  materials.  In  addition,  Zerilli 
and  Armstrong  suggest  that  a  term  may  be  necessary  in  the 
expression  to  account  for  deformation  twin  formation,  which  is 
known  to  occur  in  deformed  bcc  materials. 

3.  TWO-DIMENSIONAL  MODELING 

Once  a  suitable  constitutive  relationship  is  chosen,  and  the 
material  parameters  determined,  a  mesh,  or  grid,  model  of  the 
target  and  projectile  is  constructed.  In  this  type  of  model,  the 
mass  of  the  material  within  the  elements  is  distributed  at  the 
connecting  points  throughout  the  body. 

To  initialize  the  model,  a  velocity  is  given  to  all  of  the 
elements  in  the  model  which  describe  the  projectile.  At  some 
time  after  impact,  the  computer  program  will  calculate  the 
deformation  behavior  of  the  material  for  the  particular  impact 
velocity  chosen.  This  process  can  be  run  iteratively  until  the 
projectile  is  brought  to  rest. 

The  major  advantage  of  this  type  of  modeling  is  the  ability 
to  generate  a  historical  record  of  the  deformation  as  it 
proceeds  to  the  final  ballistic  geometry.  This  capability  of 
multidimensional  modeling  provides  a  tool  to  study  the  details 
of  the  deformation  event.  One  measure  of  success  for  the 
constitutive  relationship  is  how  well  its  prediction  of  the 
final  geometry  agrees  with  the  experimental  results. 

The  inherent  difficulties  associated  with  this  type  of 
modeling  is  the  amount  of  testing  required  to  suitably 
characterize  a  material  and  the  computational  resources 
necessary  to  run  the  model.  With  regard  to  the  former,  an 
extensive  test  matrix  must  be  constructed  and  executed  to 
adequately  characterize  the  material.  If  the  processing  history 
of  the  material  is  altered,  then  the  characterization  is  no 
longer  valid  because  of  changes  in  the  microstructure  of  the 
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material.  As  an  example,  the  parameters  determined  for  pure 
copper  tested  in  the  annealed  condition  would  not  be  equal  to 
those  found  for  copper  tested  in  a  half-hard  condition. 

In  addition,  the  constants  used  to  characterize  the  material 
are  generally  not  determined  under  the  conditions  found  in  an 
impact  event.  Therefore,  using  the  constitutive  relationships 
in  impact  studies  represents  an  extrapolation  of  the  material 
model  into  a  regime  where  little  or  no  prior  information  exists. 
Data  used  to  characterize  a  material  are  often  obtained  at 
strain  rates  several  orders  of  magnitude  below  those  found  in 
the  actual  event.  At  sufficiently  high  rates  of  deformation 
certain  physical  properties,  i.e.,  pressure,  temperature, 
density,  etc. ,  become  increasingly  important  in  the  behavior  of 
the  material.  The  temperature  terms  found  in  Equations  (41)  and 
(42)  are  attempts  to  account  for  the  thermodynamic  influences  on 
the  material  behavior. 

The  requirements  for  computer  resources  are  large  in 
multidimensional  models.  To  handle  the  extensive  amount  of  data 
which  represents  the  current  state,  or  condition,  of  the 
specimen  and  target  requires  large  core  memorv  storage  In 
addition,  thousands  of  calculations  are  necessary  to  accurately 
increment  the  process  a  small  step  in  time.  In  combination, 
these  two  basic  requirements  of  storage  and  computational  speed 
eliminate  the  personal  computer  as  an  effective  tool.  Cost  to 
the  user,  because  of  the  computer  resource  requirements,  is  the 
primary  drawback  of  multidimensional  modeling  techniques.  If 
the  factor  of  cost  can  be  overcome  then  the  use  of 
multidimensional  modeling  is  appropriate  for  the  investigation 
of  the  details  surrounding  the  deformation  process. 

In  general,  the  interest  in  multidimensional  modeling  has 
developed  simultaneously  with  that  of  one-dimensional  models.  In 
both  cases,  the  investigators  are  attempting  to  provide  materials 
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characterization  so  that  predictive  models  in  the  area  of 
target/penetrator  interaction  can  be  developed. 

Multidimensional  models  are  useful  in  target/penetrator 
interaction  for  studying  the  details  of  the  deformation  process. 
The  assumption  in  multidimensional  modeling  is  that  if  the  final 
configuration  predicted  by  the  model  is  matched  by  experimental 
results,  then  the  constitutive  relationship  is  a  success  at 
modeling  the  process.  Unfortunately,  success  at  predicting  the 
outcome  is  not  a  sufficient  criteria  for  judging  our  knowledge  of 
the  mechanisms  at  work  in  the  penetration  process.  The  transient 
nature  of  the  penetration  process  makes  it  difficult  to  obtain 
useful  data  to  confirm  the  results  of  multidimensional  models. 

For  one-dimensional  analyses,  the  measure  of  success  is  the 
ability  to  predict  the  depth  of  penetration  into  a  target 
material.  Such  analyses  can  be  performed  after  the  target  and 
penetrator  materials  have  been  characterized  by  the  Taylor  test 
method.  This  provides  the  armor/armor  penetrator  designer  with  a 
straightforward  technique  for  predicting  the  results  of  the 
target/penetrator  interaction  to  a  first  approximation  level.  By 
using  a  one-dimensional  technique,  the  effectiveness  of  a 
material  as  a  car  lidate  for  a  penetrator  or  armor  can  be 
estimated  prior  to  the  effort  and  expense  of  experimental 
verification. 
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SECTION  VII 


CONCLUSIONS 


1 .  INTRODUCTION 

This  investigation  has  produced  results  in  the  following 
areas:  evaluation  of  the  experimental  apparatus  and  procedures, 

further  insight  into  the  physical  process  resulting  from  impact, 
e--aluation  of  the  ct/0  model.  The  conclusions  resulting  from  the 
analysis  of  these  areas  can  be  described  separately,  but  it  must 
be  kept  in  mind  that  they  represent  areas  of  knowledge  developed 
in  parallel  during  the  investigation.  Improvements  to  the  test 
apparatus  and  methodology  resulted  in  further  understanding  of 
the  physical  process;  which  in  turn,  expands  the  basis  for 
evaluation  of  the  engineering  model.  Overall  the  results  show 
that  the  investigation  has  been  effective  in  improving  the 
experimental  technique  and  understanding  of  the  physical  process, 
which  can  be  applied  towards  the  enhancement  of  predictive 
models . 

2.  EXPERIMENTAL  APPARATUS  AND  METHODOLOGY 

The  engineering  design  has  succeeded  in  meeting  the  baseline 
gcal  for  producing  an  efficient  and  flexible  test  apparatus 
having  good  repeatability.  The  efficiency  of  the  design  provides 
a  cost-effective  source  for  material  characterization  data.  In 
addition,  the  test  apparatus  has  designed-in  flexibility  by 
comparison  with  systems  previously  used.  Experimentally  the  most 
important  goal  is  repeatability.  This  factor  allows  the 
structure  of  experiments  to  focus  on  extracting  information  from 
each  test,  as  opposed,  to  the  concern  that  the  test  methodology 
will  obscure  trends  in  the  data. 

Improvements  in  the  technologies  for  measuring  projectile 
velocity  and  ballistic  geometry  have  revealed  certain  trends  in 
the  data.  With  these  capabilities  the  physical  process  is  more 
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readily  understood,  and  the  possible  enhancement  of  analytical 
models  is  likely. 

3.  PHYSICAL  PROCESS  RESULTING  FROM  IMPACT 

During  this  investigation  several  features  of  the  physical 
process  were  revealed,  which  had  previously  not  been  reported  in 
the  regime  of  velocities  studied.  By  the  observation  of 
recrystallization  and  deformation  twins  in  pure  copper  specimens 
it  was  evident  that  the  initial  pressure  was  significantly  higher 
than  previously  thought.  A  model  constructed  on  the  basis  of 
shock  wave  theory  was  in  good  agreement  with  the  evidence  found 
in  the  microstructure  of  the  specimen. 

The  magnitude  of  the  pressure  found  in  materials  used  in  a 
Taylor  test  was  significant  in  that  the  basic  analysis  of  the 
process,  including  the  current  model,  uses  traditional  plastic 
and  elastic  wave  motion  as  its  foundation.  In  actuality,  there 
is  sufficient  evidence  to  indicate  that  the  deformation  in  the 
specimen  is  initiated  as  a  shock  process,  though  it  is  relatively 
short  lived  by  comparison  to  the  total  length  of  time  of  the 
event.  Evidence  which  supports  the  conclusion  that  the 
deformation  process  cannot  be  described  strictly  as  a 
plastic  and  elastic  wave  process  is  found  in  the  results  from 
the  high  speed  photography.  Data  shown  in  Figure  20  indicates 
that  to  a  reasonable  approximation  the  rigid  rod  velocity  did 
not  change  until  the  event  was  half  over.  Plastic/elastic  wave 
analysis  suggests  that  the  undeformed  end  of  the  specimen  would 
have  a  change  in  velocity  within  the  first  30  u s  of  the  event. 
This  suggests  that  the  shock  process  is  altering  the  time  frame 
of  events,  such  as  the  completion  of  mushroom  growth  and  the 
change  in  velocity  of  the  rigid  rod,  relative  to  theory  on 
plastic  and  elastic  wave  motion  in  materials. 
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The  questions  that  arise  crom  these  conclusions  are:  when 
does  the  shock  wave  attenuate  into  plastic  and  elastic  waves; 
what  is  the  shape  of  the  wave  front,  in  the  specimen,  which  is 
considered  planar  for  simplicity  in  the  analysis?  The  first  of 
these  questions  could  be  addressed  by  examining  the  rate  at 
which  a  shock  wave  is  attenuated  as  it  moves  through  a  material. 
In  addition  to  attenuation  by  the  motion  of  the  shock  front,  the 
free  boundary  of  the  material  produces  a  release  wave  which  will 
have  some  influence  on  the  process.  The  second  question  could 
be  addressed  by  a  careful  examination  of  the  microstructure 
resulting  from  impact.  Since  shock  waves  are  produced  by 
pressure  levels  that  have  certain  effects  on  the  microstructure, 
it  may  be  possible  to  determine  the  shape  and  extent  of  the 
shock  process  by  a  thorough  microstructural  examination. 

4.  a/£  MODEL 

The  Taylor  test  was  developed  as  a  simple  method  of 
obtaining  data  to  characterize  the  strength  of  a  material  under 
impact  conditions.  The  a/p>  model,  Reference  8,  was  successful 
as  an  analysis  which  accounted  for  mass  transfer  from  the  rigid 
portion  of  the  rod.  This  analysis  was  constructed  based  on  the 
assumptions  that  the  constitutive  relationship  was  perfectly 
plastic,  that  radial  inertia  effects  could  be  ignored,  and  that 
the  process  is  described  by  plastic  and  elastic  wave  motion. 

The  assumption,  regarding  the  constitutive  relationship  and  the 
effects  of  inertia,  characterize  the  analysis  as  a  first 
approximation  suitable  as  an  engineering  model. 

The  term  engineering  model,  as  used  here,  is  defined  as  an 
analysis  which  is  based  on  the  assumption  that  an 

interrelationship  exists  among  the  test  variables,  i.e.,  density, 
velocity,  geometry,  etc. ,  that  will  reveal  the  strength  of  a 
material  on  impact.  The  fact  that  the  analysis  was  developed  as 
an  engineering  model  must  be  considered  when  evaluating  its 
success . 
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A  recently  proposed  model  of  high  strain  rate  behavior  in 
pure  copper  provides  a  measure  of  success  for  the  current 
investigation.  A  stress/strain  curve,  Figure  36,  generated  by  a 
model  proposed  by  Follansbee  and  Kocks ,  Reference  17,  shows  the 
predicted  behavior  of  pure  copper  under  high  strain  rate 
conditions.  Given  the  level  of  approximation  found  in  the  *//5 
model,  it  has  provided  data  which  is  in  good  agreement  with  the 
strength  of  pure  copper  given  by  the  model  of  Follansbee  and 
Kocks.  At  the  level  of  approximation  for  which  it  was  designed, 
the  cc/0  model  is  successful  at  predicting  the  strength  of  a 
material  under  impact  conditions. 


Figure  36.  Stress/Strain  Response  of  Pure  Copper 
Under  High  Strain-Rate  Conditions 

The  information  obtained  from  the  Taylor  test  and  the 
analysis  developed  by  investigators  Jones,  Gillis,  and  Foster, 
Reference  8,  will  provide  the  armor/armor  penetrator  designer 
with  a  relatively  simple,  yet  effective,  tool  for  the  estimation 
of  raterial  strength  at  high  rates  of  deformation.  While  the 
model  is  not  intended  to  describe  the  detailed  behavior  of  the 
physical  process,  it  does  provide  a  first  approximation  approach 
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to  determining  the  material  property  necessary  in  the 
development  of  terminal  ballistic  models. 
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APPENDIX  A 
TEST  PROCEDURES 
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APPENDIX  A 


TEST  PROCEDURES 


1)  Remove  target  housing  cover.  Clean  debris,  if  any, 
off  of  the  plexiglass  windows. 

2)  Remove  and  clean  pressure  transducers,  if  necessary. 

3)  Prior  to  mounting  the  pressure  transducers  in  the 
muzzle,  apply  a  thin  layer  of  silicon  grease  over  the 
piezo-crystal  surface  to  act  as  a  damper  on  the  pressure 
wave . 

4)  Check  the  voltage  output  from  the  infrared  detector 
circuit.  A  low  voltage  state  indicates  debris  is 
covering  the  detector  lens  and  must  be  removed. 

5)  Clean  the  infrared  detector  lens,  if  necessary. 

6)  Rotate  the  target,  if  necessary,  and  position  the 
fiducial  (magnet)  on  the  target  face  one-half  inch  above 
the  point  of  impact. 

7)  Bore  sight  the  high  speed  camera,  or  the  framing 
camera,  if  necessary.  The  camera  lens  should  be  as  near 
parallel  to  the  target  face  as  possible. 

8)  If  high  speed  photography  is  to  be  used,  load  the 
camera  with  film. 

9)  Set  the  desired  standoff  distance  between  the  muzzle 
and  the  anvil.  Tighten  the  bolts  on  the  v-block  mounts 
to  secure  the  launch  tube. 

10)  Replace  the  target  housing  cover. 

11)  Weigh  specimen,  (grains) 

12)  Measure  specimen  diameter,  (inches) 

13)  Weigh  propellant  charge,  (grains) 

14)  Load  the  measured  quantity  of  propellant  into  a 
cartridge. 

15)  Insert  a  ball  of  cotton  into  the  cartridge  case  and 
pack  against  the  propellant. 

16)  Obtain  the  safe/arm  panel  keys. 

17)  Insert  a  specimen,  followed  by  an  obturator,  into 
the  breech  end  of  the  launch  tube. 
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18)  Position  the  specimen  and  obturator  at  the  correct 
depth  in  the  launch  tube  using  depth  measurement  gauge. 

19)  Insert  the  cartridge  case  behind  the  specimen  and 
obturator.  Screw  the  breech  cap  on  the  barrel  until  the 
scribed  marks  are  aligned. 

20)  Attach  the  firing  solenoid  to  the  breech. 

21)  Remove  shunting  connector  from  the  firing  line  and 
connect  the  firing  line  to  the  firing  solenoid. 

22)  Insert  the  firing  pin  to  the  prescribed  depth  into 
the  firing  solenoid. 

23)  Remove  shorting  plug  from  firing  line  and  connect 
firing  line  to  firing  power  supply  output  panel. 

24)  Arm  the  firing  power  supply  circuit. 

25)  Arm  the  firing  signal  control  panel. 

26)  After  arming  the  firinq  signal  control  panel,  moni¬ 
tor  the  cnarge  level  of  the  firing  power  supply  until 
the  predetermined  voltage  level  is  reached.  Once  the 
appropriate  energy  level  is  reached  the  firing  count¬ 
down  can  begin. 

27)  At  t  minus  20  seconds,  the  multitrack  signal 
recorder  is  turned  on. 

28)  At  t  minus  5  seconds,  when  using  the  high  speed 
movie  camera,  the  photographic  light  is  turned  on. 

29)  At  t  minus  3  seconds,  the  high  speed  camera  is 
turned  on. 

30)  At  t  equal  to  0,  the  firing  signal  is  sent  to  the 
high  speed  movie  camera  and  to  the  firing  solenoid. 

31)  Disarm  the  control  panel. 

32)  Disarm  the  firing  circuit. 

33)  Remove  the  firing  circuit  wire  connected  to  the 
firing  solenoid. 


34)  Unscrew  the  breech  end  cap  and  extract  the  spent 
cartridge  case. 

35)  Clean  debris  from  the  launch  tube  bore  with  a  gun 
cleaning  rod  tool. 
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36)  Remove  target  cover  and  retrieve  the  impact 
specimen  and  obturator.  Discard  obturator. 

37)  Remove  high  speed  movie  film,  or  framing  camera 
film,  if  necessary. 

38)  Measured  from  the  oscilloscope,  record  the  duration 
of  time  between  the  leading  edge  of  the  pressure  trans¬ 
ducer  signals. 

39)  Obtain  a  strip  chart  output  of  the  transducer  sig¬ 
nals  from  the  multichannel  recorder. 

40)  Using  the  strip  chart  output,  digitize  the  reference 
signal  (frequency  =  100  khz).  Digitize  the  leading  edge 
of  the  pressure  transducer  signals. 

41)  Obtain  a  print  out  of  the  digitized  data  converted 
to  specimen  velocity. 

42)  Measure  the  final  diameter  of  the  mushroomed  portion 
of  the  specimen. 

43)  Measure  the  final  length  of  the  specimen. 

44)  Measure  the  undeformed  length  of  the  specimen. 
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MATERIAL  OFHC  COPPER 
YIELD  STRENCTH:  M 
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MATERIAL:  OFHC  COPPER  (Com. 
YIELD  STRENGTH:  Mpa 

ULIMATE  STRENGTH:  304  Mpa 
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MATERIAL:  DPTE  COPPER 
YIELD  STRENGTH:  Mpa 

ULTIHATE  STRENGTH:  300  Mpa 
DENSITY:  8.94  CRAMS/CUBIC  CENTIMETER 
PILE  WAS  UPDATED  ON  DECEMBER  16,  1988 


a)  EXCESSIVE  DEEGRMATiON 

b)  OBLIQUE  IMPACT 


MATERIAL:  6061  T6  ALUMINUM 
YIELD  STRENGTH:  315  Mpa 
ULTIMATE  STRENGTH:  340  Mpa 
DENSITY:  2  71  CRAMS/CUBIC  CENTIMETER 


MO 

Co 

MO 

06 

CM 

ao 

mo 

O' 

vo 

Cm 

06 

ao 

s© 

— c 

*c 

in 

CM 

n- 

CD 

06 

Os 

CO 

as 

VO 

rs 

>a 

-■A 

CO 

(2 

© 

r-x 

MO 

M3 

© 

C6 

M) 

ao 

o 

a 

n 

>a 

I/O 

in 

ao 

Ul 

Ul 

U4 

CM 

O' 

CM 

<y 

GO 

co 

CM 

CM 

CM 

© 

•a 

O' 

O' 

© 

H 

CM 

>o 

© 

n 

O 

<S| 

MO 

MT 

ao 

o 

M3 

A-4 

rv 

CM 

MJ 

r— 1 

in 

CO 

co 

© 

O 

cC 

(S 

sO 

M3 

CO 

ao 

CO 

MO 

•a 

mo 

I/O 

-a 

-a 

CO 

n 

O' 

<r 

in 

O' 

<r 

a- 

-a 

•a 

ST 

xT 

■a 

1  CM 

CM 

CM 

Me 

M 

MO 

O' 

ao 

r*. 

CM 

Me 

04 

<r 

ao 

r*». 

rx 

CM 

ao 

M3 

Ml 

CO 

CO 

rx. 

cs 

CM 

CM 

CO 

CM 

xj 

Ml 

CO 

ao 

M3 

ao 

CM 

I/O 

Mt 

© 

ec 

co 

ao 

mo 

ao 

CM 

O 

Mt 

•a 

O 

I/O 

M3 

CO 

CO 

CO 

'O* 

in 

Mt 

O' 

ao 

O' 

M3 

© 

o 

Mi 

f*x 

—> 

O' 

as 

mo 

<r 

r- i 

CO 

Ml 

— X 

o 

i/O 

i/O 

MO 

MO 

■M 

o 

Me 

•-4 

M3 

CM 

O' 

oO 

1—4 

Ml 

o 

sO 

O 

as 

o 

O 

Ml 

© 

O' 

ao 

o 

rx. 

CM 

Me 

Ml 

4-4 

*”4 

Ml 

t-4 

r*. 

CM 

* 

^4 

CM 

CM 

rM 

*-4 

«— 4 

Ml 

CM 

l>4 

CM 

4-4 

CM 

CM 

* 

CM 

Mt 

Mt 

Mi 

Me 

Me 

Ml 

M3 

-X. 

O' 

Mt 

CO 

MO 

co 

— 

CO 

Me 

SO 

-x 

vO 

fm4 

o 

Os 

rs 

CM 

Ml 

CO 

as 

■a 

CO 

as 

MO 

•— • 

mo 

CM 

Mt 

n 

co 

co 

ao 

CO 

—4 

Mt 

ao 

SO 

CO 

CM 

co 

n 

O' 

vO 

X 

co 

~x 

M3 

CM 

O 

M3 

CM 

rx 

Me 

co 

r-» 

^4 

O' 

CM 

Mt 

C 

•x. 

~x 

CM 

CM 

rs 

O' 

ao 

Mt 

ao 

CM 

AX. 

I/O 

© 

CM 

CM 

Me 

M3 

CM 

CM 

CO 

CO 

M3 

Mt 

CO 

•a 

CM 

O 

1-4 

CM 

■a 

MC 

CM 

<*• ' 

CO 

a 

CM 

CO 

CO 

M« 

Ml 

Mt 

© 

O 

O 

Ml 

© 

© 

© 

© 

© 

© 

Mt 

C 

© 

© 

© 

© 

C 

© 

© 

Mt 

© 

o 

© 

o 

o 

T-) 

It 

Me 

Mt 

Me 

Mt 

Ml 

ao 

CO 

CM 

_4 

ao 

rx. 

M3 

Cx 

Os 

_4 

CM 

CO 

03 

sO 

o 

<T 

I/O 

M3 

CM 

CO 

M3 

r»x 

Ax. 

CO 

00 

ao 

<T 

O' 

CO 

M3 

-a 

Os 

CM 

CO 

O' 

GO 

n 

CM 

mO 

n 

xT 

CO 

© 

as 

© 

ao 

O 

CM 

xT 

MT 

as 

M3 

CO 

in 

co 

O 

V? 

o 

mo 

Px. 

rx. 

i/O 

xj 

n 

px. 

CO 

CM 

CM 

CO 

xT 

<T 

CD 

-T 

<r 

mo 

M3 

© 

•a 

O' 

r—* 

© 

© 

i/0 

r*x 

xt 

1—4 

© 

.— 1 

CO 

03 

Cv 

CM 

M3 

1 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

A*0 

ooQOOQaoao,s^r,^r^*"^r^-'OsOvC'/>i/5j^vy^^roOOOO’,^~^r'>' 

4<T«J<T'J!>0'l7'0'05<ra'C5^0'J»^0\<J>^vr^sT>*C0CDaD 


J  U1  c 
u J  o  — 

X  QC  ^3 

3  <  »- 
oc  z  © 

Om  CJ 


.  CM  CM  CM  iT  CO  CQ  i*1  co  n  O'  O'  MO  <— »  — »  -x  MO  CO 
^M-*NNN-'»>*0«MNrY^CN 


M  C5  HA  !*» 


X 

H 

o  CP 

q  z  O- 

J  U]  X 
w  ac 
>— *  H 
>-  </> 

(X 

U) 

H  v> 
J  Ul  u 
<  X  Z 
z  3  o 
—  ~  z 
u.  o  »« 


<3  r-4  O'  + 

CM  O'  Ml  *-  ^  •(jNHNOv.t'^'ONrvN^ 

•  *  n  it  «  -mo . Me 

in  n  ao  M  O^'JClO'fl'fl'OOrvCsjMCOM 

n»  so  ao  Mt  vjCi>rm^c^«jcotN'omo^M 
^  co  M  Co  co  co  c>  co  co  co  co  co  co  co  Me 

*  * 

*  * 

to  ^  as  ■*  ^tN0)fovfl(N'0^'0Mv0in(00>^M 

•  •  * . * 

r-  ro  cm  Me  iflO'0^*J>f^«CN^>?rorsa0CSM 
rs  in  »  *  HO^'-O^c'Ofonoroifl'OHit 
>o  vfl  in  ♦  in'Dmininyjcoinnn^in^nin* 
Ml  Mi 

Mt  Me 

Ml  Ml 

Mi  ao  O'  *  cn  n  c  <?>  O'  o.  •*  otvntTiciHsjm* 

Ml  ao  03  Ml  o  CO  O  H  «J  h  Mt  C'0OM>}O'CO'JiAMl 

Mt  CO  CO  Ml  in  (O  4  c  n  Mt  ^  «T  <r  >J  >J  O  -Ml 

Ml  -Ml . Ml . OOMt 

MtOOMiOOOOOOMiOCOOOO  Me 


CM  xj  a  M*  C3''0'OsJNCOvono-^)cO(\inCN* 

JO  o  M  OiH'JCNCNeOO'CNlT-'OOCNT'COMI 

CM  CM  -4  Ml  CtNtNWtSi-CinMnvJ^Oin^* 


o  f  O  *  OOOOOOOO 


O  ©  O  O  O  O  Mi 


(n  Ch  in  m  oo  o  sj 

CO  fM  1-0  <—»  fx.  ^ 

iO  in  n  ®  h  n  n 

oi  in  co  in  n  >o  o 

(M  po  o  n  ^  in  n 


in  v©  so  vo  co  mo  cm 

O  to  O  H  N  N  (Nl 

n  h  n  cn  rs  m  n 

n  n  n  n  n  n  n 


•J  ^  C  Jl  Jl  \  (\ 

>J  Y3  «J  sj  CN  lO  O 

^  ^  n 

•  o 

O  O  C  ©  o  c 


n  n  o  -D  ®  (”  m 

O'  lO  CO  ®  M  ^ 

•«  00  is  in  -  — 

O  •  ■  ■  ■  —  ■ 

o  ©  ©  — «  — 


n  •«  >4  «  vOiono''00''OinsnO'P^O'>>.0'M 
(C  o  O'  Ml  (Mr'NOO'0(N'0(N-J-«0''i)®s< 
^cciMtnininino«jaoNooO  ■  c>  ?m  — <  Mt 

■Ml . ©•••<** 

OOOM  OOOOOCOOC--0  mi 


ao  ao  ao  co  «  m  cn 

>y  O  O  cm  ^  mo  ao 

O'  O'  o>  ao  o  ■  -~r 


NaocN'Ot,n\otnsiniocino*j(\"Jrsco'OP5Hn'0^r'oO' 

aooaofniOin'flntnvOflO'flnocYNfsi-NP'ioinOinfsrso 

ONjO'aOHfn^inNorscO'fictNO'oninncicnjnH 

NNMflC'HNrt'fltvN^sNwaO'O^O'intn^inNinins 

HN<MW^CMNCNN(n»R<NNNNniHtNN  CN  CM  CM  CM  0*4  CN  CN 

^(n-oncwc.vo^in^in^vO'6von‘«-N.O'C'Mrs,^^^^ 
■  ■  CO  >0  O  ■  •— *  CM  V/O  ■  rx.  ■  O 

ifirs©inc^(ns«oOCNO-aoonc«inrsn(noio®oi 

O'  cx  ao  O'  co  h  <^)  ci  »j  n  O'©--*  O'  cs  \  n  co 

CO  N  N  O  to  I/O  ao  N  rs  33  O'  n  r«x  ao  rs  ao  ao 


l/OMOMOMO«/OCMCMfMCM<NCOOCO<f>y.aMOi/OmMOMOMOnin©0©  I 

f*  ~4  ~  r*  s  N  N  N  s  ^  H  j 

4a001'0'OOc^N«n®0'0'tNrsO'srs«<ycin^f^^®Oi  I 

CMtYvO'ONnniicH^<r-*i^'J(N<Hio-4yj^^^(\<is.<fl„  { 


96 


a)  EXCESSIVE  DEFORMATION 

b)  OBLIQUE  IMPACT 

c)  PRIMER  FAILURE 


MATERIAL  2024  T4  ALUMINUM 
YIELD  STRENGTH:  400  Mpa 
ULTIMATE  STRENCTH :  500  Mpa 

DENSITY:  2.78  GRAMS/CUBIC  CENTIMETER 
FILE  WAS  UPDATED  DECEMBER  16,  1988 
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Figure  D-l  .  Load  Response  of  0F11C  Copper 
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Figure  U-2 .  Load  Response  o£  DR'J’L  Copper 


INCREASING  TIME 

LOAD 


Figure  D-4 .  Load  Response  of  2024-T4  Aluminum 


APPENDIX  E 


Figure  E-l.  Schematic  Diagram  Showing  the  Parameters 
Used  in  the  Hugoniot  Model 

Nomenclature : 

P  =  Pressure  at  the  target/specimen  interface 
PI  =  Initial  density  of  the  target 
p-  =  Initial  density  of  the  specimen 

p,  =  Density  of  the  target  material  behind  the  shock 
front 

p 2  =  Density  of  the  specimen  material  behind  the  shock 

front 

u  =  Target/Projectile  interface  velocity 
Uis  =  Shock  wave  velocity  in  the  target 
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ulp  =  Particle  velocity  behind  the  shock  wave 
in  the  projectile 

u2s  =  Shock  wave  velocity  in  the  specimen 
u2p  =  Particle  velocity  behind  the  shock  wave 
*  in  the  specimen 
v  =  Impact  velocity 

L  =  Original  length  of  the  specimen 

A  =  Cross-sectional  area,  assumed  constant 

At  =  Time  increment 

To  predict  if  deformation  twins  can  form  upon  impact  of  pure 
copper  test  specimens,  a  model  was  developed  to  predict  pressure 
at  impact.  The  pressure  can  be  determined  from  a  set  of  equa¬ 
tions  derived  from  impulse-momentum  and  continuity  equations, 
together  with  Hugoniot  relationships.  It  is  assumed  for  this 
analysis  that  the  velocity  of  impact  is  sufficient  to  generate 
shock  waves  in  both  the  target  and  the  projectile.  The  analysis 
is  that  used  in  the  classical  flyer  plate  experiment,  but  it  is 
only  valid  at  impact,  or  for  a  short  time  afterwards,  in  the 
Taylor  test. 

At  the  free  boundary  of  the  specimen,  radial  motion  of  the 
material  undergoing  deformation  will  invalidate  the  conditions  on 
the  model.  Therefore,  the  model  represents  only  the  initial 
conditions  at  impact.  However,  near  the  axis  of  the  specimen  a 
state  of  high  pressure  exists  until  radial  release  waves  reach 
this  region  from  the  free  surface.  This  region  of  the  specimen 
was  the  primary  focal  point  of  the  microstructural  investigation. 

Conservation  of  mass  can  be  used  to  establish  a  relationship 
between  the  various  velocity  terms  and  the  density  of  the 
material.  For  the  target,  the  mass  under  consideration  is  that 
volume  of  material  that  decreases  because  of  the  shock  front. 
Figure  E-l.  The  original  mass  of  material  can  be  equated  to  the 
final  mass,  at  some  incremental  time  after  impact,  as  follows: 

/5lulsAtA  =  ?iA(Uis  -  ulp)At 
Dividing  both  sides  by  the  common  terms  gives 


(E-l) 


Pluls  -  Pl<Uis  ~  ulp) 

The  particle  velocity,  u^p,  in  the  target  is  equal  to  the 
interface  velocity,  u.  After  substituting  for  the  particle 
velocity,  Equation  (E-l)  can  be  written  as 

PlUis  =  pi(Uis  ”  u) 

This  expression  states  that  the  mass  of  material  must  be 
conserved  before  and  after  a  shock  wave  moves  through  a  volume 
of  material. 

For  the  projectile,  a  similar  relationship  can  be  expressed 
as  follows: 

P2U2s  =  P2(u2s  +  u2p)  (E-2) 

In  the  projectile,  the  particle  velocity,  u2p,  is  related  to  the 
impact  velocity,  v,  and  the  interface  velocity,  u,  by  the 
expression 

u2p  =  v  -  u 

Substituting  for  u2p,  Equation  (E-2)  becomes 
P2u2s  =  P2(u2s  +  v  -  u) 

Using  the  impulse  equation  for  the  same  volume  of  material 
in  the  target,  a  relationship  is  derived  between  the  velocity  of 
the  material  and  the  pressure: 

PAAt  =  p1A(Uls  -  U)UAt  (E-3) 

Dividing  both  sides  by  the  common  terms  gives 

p  =  Pl(uls  -  u)u 

This  expression  describes  the  relationship  between  the  velocity 
variables  and  the  pressure,  P.  For  the  projectile,  the 
relationship  can  be  expressed  as 

P  =  ~P2  (^2s  +  v  "  u)u  “  P2u2sv  (E-4) 
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Equations  (E-l)  through  ( E— 4 )  contain  six  unknown 
parameters:  P,  u,  U]_s,  Tj 2 s,  ~pi,  and  p2  •  From  shock  experiments, 
two  additional  expressions  can  be  obtained,  giving  a 
relationship  between  the  shock  velocities,  Uls,  and  U2s,  and  the 
particle  velocities,  ulp,  and  u2p.  The  data  from  these 
experiments  can  be  fit  to  a  straight  line  approximation  of  the 
type 


Uis  -  +  biuip  -  aj_  +  b^u 


(E-5) 


where  a^  and  b^  are  quantities  determined  from  shock  experiments 
of  the  target  material,  Reference  E-l.  A  similar  expression  can 
be  used  for  the  specimen  material: 


l2  +  b2u2p  =  a2  +  b2  (v  -  u) 


(E-6; 


Solving  for  the  pressure  requires  that  the  interface  velo¬ 
city  u,  be  determined.  The  unknown  densities,  plf  and  p2 ,  must 
be  eliminated  in  Equations  (E-3)  and  (E-4)  by  substituting  from 
Equations  (E-l)  and  (E-2) .  The  shock  velocities  can  also  be 
eliminated  in  Equations  (E-3)  and  (E-4)  by  substituting  from 
Equations  (E-5)  and  (E-6) .  The  pressure  equations  can  be  set 
equal  to  one  another  and  algebraically  manipulated  into  Lhe 
following  form: 


where 


+  Mu  +  N  =  0 


—  (73]_  +  a2  +  2b2v) 


b2  -  7^1 
(a2  +  b2v)v 
b2  _  7bi 


so  that 


T  =  P\/P2 


-M  ±  lMz  -  4  (N) 
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Substituting  the  interface  velocity,  u,  into  the  equations 
the  unknown  quantities,  Uls,  U2s,  pi,  ~p2>  and  P'  can  be 
obtained. 

A  sample  calculation  using  a  steel  4340  target  material  with 
a  copper  projectile  moving  at  200  m/s  predicts  an  interface 
velocity  of  100  m/s.  This  value  of  u  can  be  substituted  back 
into  the  equation  set  to  give  a  pressure  value  of  3.9  GPa.  This 
pressure  level  is  almost  an  order  of  magnitude  higher  than  the 
high  strain  rate  flow  stress  of  the  material,  which  is  about  0.5 
GPa.  According  to  Johari  and  Thomas,  Reference  E-2 ,  pressure 
greater  than  2  Gpa  is  sufficient  to  cause  the  formation  of  twins 
in  copper  material. 
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APPENDIX  F 


RAW  DATA  SHEETS 

TAYLOR  IMPACT  TESTS 

TEST 

NO _ 

OATE _ 

LAUNCHER  DESIGNATION _ 

STANDOFF  DISTANCE _ 

TARGET  DESCRIPTION:  MATERIAL _ 

0 IMENS IONS _ 

HARDNESS _ 

OBLIQUITY _ 

SPECIMEN  DESCRIPTION:  MATERIAL _ SOURCE 

LENGTH _ 

DIAMETER _ 

WE  I GHT _ 

OBTURATION  WEIGHT _ 

TOTAL  WEIGHT _ 

PROPELLANT  _ WEIGHT _ 

TEMPERATURE _ 

HUMIDITY _ 

VELOCITY _ 

CAMERA _ 

TYPE  FILM 


FPS _ 

F/STOP 
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